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2 ABSTRACT 

-4— > 

The recent discovery of the "weak field, old magnetar", the soft gamma repeater 
> ■ SGR 0418+5729 , whose dipole magnetic field, Bdip, is less than 7.5 x 10^^ G , 

oo! 

Q>^ , has raised perplexing questions: How can the neutron star produce SGR-like bursts 

. with such a low magnetic field? What powers the observed X-ray emission when 

neither the rotational energy nor the magnetic dipole energy are sufficient? These 
observations, that suggest either a much larger energy reservoir or a much younger 
true age (or both), have renewed the interest in the evolutionary sequence of mag- 
netars. We examine, here, a phenomenological model for the magnetic field decay: 
Bd\p oc B]^^ and compare its predictions with the observed period, P, the pe- 
riod derivative, P, and the X-ray luminosity, Lx, of magnetar candidates. We find a 
strong evidence for a dipole field decay on a timescale of ^ 10'^ yr for the strongest 
(Sdip ^ 10^^ G) field objects, with a decay index within the range 1 ^ a < 2 
and more likely within 1.5 ^ a < 1.8. The decaying field implies a younger age 
than what is implied by P/2P. Surprisingly, even with the younger age, the energy 
released in the dipole field decay is insufficient to power the X-ray emission, suggest- 
ing the existence of a stronger internal field, Bi^t . Examining several models for the 
internal magnetic field decay we find that it must have a very large (> 10^^ G) initial 
value. Our findings suggest two clear distinct evolutionary tracks - the SGR/AXP 
branch and the transient branch, with a possible third branch involving high-field 
radio pulsars that age into low luminosity X-ray dim isolated neutron stars. 
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1 INTRODUCTION 



Soft Gamma Repeaters (SGRs) and Anomalous X-ray Pulsars (AXPs) are two classes of pulsating 
X-ray sources thought to be highly magnetized Neutron Stars (NSs) whose high-energy emission 
is powered by the dissipation of their magnetic field. Hence they are deemed magnetars. They 
have large rotation periods (5 s < P < 12 s) with relatively large time derivatives (P ~ 10^^^ — 
10^^° s s^^), implying (through magnetic dipole breaking) large surface dipole field strengths 
(Pdip ^ 3 X 10^'^ G or energy -^Bdip ~ 2 x 10^^ erg) and relatively young spindown ages (Tc ~ 
10^ — 10^ yr). Thus, the decay of their dipole fields on the timescale of their spindown ages 
appears to be capable of accounting for their large persistent X-ray luminosities, of Lx ~ 10^^ 
erg s~^ ([Thompson & Duncanlll995|) . On the other hand, rotational energy losses are 1-2 orders 
of magnitude too low to account for their measured Lx, and the absence of binary companions 
strongly argues against accretion as a power source. 

Both SGRs and AXPs can emit sporadic, sub-second (~ 0.1 s) bursts of hard X-rays to soft 7- 
rays, which release ~ 10^^ — 10'^^ erg at typically super-Eddington luminosities. Evolving stresses 
from a decaying magnetic field > 10^"^ G can in principle stress the rigid lattice of the NS crust 
beyond its yielding po int, leading to sudden release of stored magnetic energy, which is believed 



to trigger these bursts ([Thompson & Duncan 



19961 : 



Pema & Pons 



2011 ). 



AXPs tend to be less burst-active, and have somewhat weaker dipole fields and larger spindown 
ages compared to SGRs. This l ed to the suggestion tha t SGRs evolve into AXPs as they age 



and their magnetic field decays (|Kouveliotou et al. 



19981). The effect of the decay of the dipol e 



magnetic field on the spin evolution of magnetars was previously considered (|Colpi et al 



200' 



but its ultimate implications for their expected X-ray emission have not yet been fully explored. 
In addition to classical SGRs and AX Ps, the magnetar fam ily includes also transient sources. 



First discovered in 2003 (XTE Jl 8 10- 197; 



Ibrahim et al. 



20041) . the group of transients is now the 



largest among magnetar candidates and is rapidly growing, thanks to the improved detection ca- 
pabalities of the Fermi Gamma-ray Burst Monitor (GBM). The distinctive property of transients 
is that they are discovered only thanks to outbursts, when they emit typical magnetar-like bursts 
accompanied by a very large increase (~ two orders of magnitude) of their persistent emission. 
Their timing parameters can only be measured in outburst and they match well those of persis- 
tent SGRs/AXPs. The implied rotational energy losses are much lower than the outburst-enhanced 
persistent emission. However, when the much weaker quiescent emission of transients could even- 
tually be measured, it was found to be below the level of rotational energy losses. 
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Typical magnetar-like bursts and outburst was also detected from the 0.3 s , allegedly rotation 



powered X-ray pulsar PS R J1846-0258, with a dipole field of 4.8 x lO^'^ G (iGavriil et al 



2008 



Kumar & Safi-Harb 



2008h . and the radio pulsar PSR 1622-495 0, with a spin peri od of 4.3 s and a 



dipole field of 3 x 10^"^ G, displayed a flaring radio emi ssion (iLevin et al 



very similar to those of two transient radio magnetars (|Camilo et al 



20101). with properties 



2006 



2007b. On the other 



hand, there are radio pulsars with dipole fields comparable to the weaker field magnetars, and 
larger than PSR Jl 846-0258, but showing no sign of peculiar behaviour. Thu s, a strong d ipole 



field does not seem a sufficient condition for powering magnetar-like emission (jKaspi 



2010h . 



The greatest surprise, however, came from the recent discovery of SGR 0418-1-5729 by the 
Fermi Gamma-ray Burst Monitor . This source, previously undetected in X-rays, was discovered 
on 5 June, 2009, through the emission of two distinct, sub-second magnetar-like bursts of low 



luminosity (< 10'^'^ erg s ^) and total energy ~ 2 and 4 x 10'^'^ erg, (|van der Horst et al. 



37 



2010h . 



Following the bursts it became detectable as a bright, pulsating X-ray source with a period of 
~9.1 s and a luminosit y of Lx ~ 10^*^ erg s~^ that decayed by 2 orders of magnitude during 
the following 6 months (lEsposito et al.ll2010|) . However, only an upper limit could be put on its 
period derivative, P < 6 x 10~^^ s s~^, translating to an up per limit of -Bdip < 7.5 x 10^^ G on 



its dipole field strength and a spindown age > 24 Myr (|Rea et al. 



2010 ). This demonstrates 



that magnetar-like activity can be present in NSs with rather standard dipole magnetic fields. This 
finding represents a breakthrough in our understanding of magnetars, as it was unexpected that 
bursts could be produced in such a low-field object, if they are associated with sudden crustal 
fractures. 

Moreover, it is clear that the dipole field of SGR 0418-1-5729 cannot power its X-ray emission 
if it decays on the timescale of its spindown age (of Tc > 24 Myr). The latter would require a much 
stronger field, of > 5 x 10^"^ G, in order to po\ yer the weakest level of X-ray emission at which 



this source was found, 1.5 yr after the outburst (IRea et al 



201ol). That i s, either its total magnetic 



energy is >5000 times larger than that the dipole (|Turolla et al.ll201 1|) . or its dipole field decay 
time is >5000 times smaller than its spindown age (or some combination of the two). 

Bearing on the above arguments, we address, in this work the power source of magnetar can- 
didates and their spin evolution as their dipole fields decay. The latter are jointly compared with 
the observed properties of the full sample of magnetar candidates, thus enabling us to draw robust 
conclusions. In §[2] we describe the main properties of the different classes of objects of interest 
and our observed sample. In §[3] we provide a simple analytic formalism for the spin evolution of 
NSs whose dipole field decays as Bdip = — -Bdip/^d oc -B^J", and discuss its general properties. A 
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summary of the mechanisms for magnetic field decay in NSs is provided in §|4l In §[5] we show that 
there is indeed a strong evidence for effective decay of the dipole fields of these sources, which is 
~ 10'^ yr timescale for the strongest fields of SGRs and AXPs, and provide quantitative constraints 
on the most likely decay mechanisms (we show that 1 < « < 2 is required). These constraints are 
made tighter (1.5 < a < 1.8) in § |6l where we show that decay of the dipole component alone 
is insufficient to power the X-ray emission and that a stronger power source is required, presum- 
ably a stronger internal magnetic field. We examine two models for internal magnetic field decay 
in § HI Comparing with observations we derive basic constraints on the initial values and decay 
properties of such a component. Our conclusions are summarized in § [H Finally, in § [9] we build 
on our conclusions to outline a self-consistent evolutionary scenario for magnetar candidates and 
speculate about their relation to other classes of objects. 



2 SOURCE CLASSES AND THE OBSERVATIONAL SAMPLE 



In order to verify the role of field decay in NSs with strong magnetic fields we consider here the X- 
ray and timing properties of magnetar candidates. These are usually identified by their measured 
values of P and P (hence the inferred dipole field, -Bdip, and spin-down age, Tc), detection of 
peculiar burst/outburst activity and persistent X-ray luminosities that exceed rotational energy 
losses. Sources are divided into two main groups, largely for historical reasons. SGRs are typically 
identified by the emission of trains of sporadic, sub-second bursts of X-to-7-ray radiation with 
super-Eddington luminosities (~ lO'^^ — 10'^^ erg s^^). More rarely they show much more powerful 
events, the Giant Flares (GFs). These are typically initiated by a spike of 7-ray emission lasting 
half a second and releasing ~ 10'^'' — 10"^^ erg at luminosities ~ 10"^^ — iC^ erg s^^ (assuming 



(~ 10^1 - 10^2 g-i) markedl 


y pulsed at the NS spin 


Mazets et al. 


1979; 


Cline et al. 


1998 


Hurley et al.l 


1999; 


Terasawa et al. 


2005; 


Palmer et al. 


2005 


). 



The highly super-Eddington luminosities of bursts and GFs definitely rules out accretion as a 
viable power source. Indeed, it first hinted to the key role of superstrong magnetic fieldil] which 
reduce the electron scattering cross-section much below the Thom son value. In par ticular, B > 



3 X 10^^ G were required to explain luminosities up to 10''^ erg s ^ (|Paczynski 



1992 ). The several 



minute-long pulsating tails are very similar in all three GFs detected so far. They suggest that a 



^ It is customary to adopt, just as a reference value for the strength of magnetar fields, the critical field -Bqed ~ 4.4 X 10^'' G, at which the 
energy of the first excited Landau level of electrons equals their rest-mass energy 
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remarkably standard amount of energy, ~ a few x 10^'' erg, remains trapped in a small-size, closed 
field line region in the magnetosphere (~ a few stellar radii), lik ely in the form of a pair-pho ton 



plasma, and slowly radiated away (the "trapped fireball" model of [Thompson & Duncan 



1995|). A 



magneto spheric field > a few x 10^"^ G would naturally account for this. 

Anomalous X-ray Pulsars were typically identified through the peculiar properties of their per- 



sistent emission ( Mereghetti & Stella 



19951). However, the distinction between these two classes 



has come into question since SGR-like bursting acti vity has been detected in many APXs (IGavriil et al 



20021 : 



Woods & Thompson 



20061 : 



Mereghettill2008|) . To date, AXP have never displayed a GF and 



are, in general, characterized by less frequent, and somewhat less energetic, burst activity. 

There are persistent and transient sources both in the SGR and the AXP classes. The emission 
from the persistent sources is known to be pulsed at the NS spin period and to be characterized by 
a blackbod y (BB) cornponen t, with kT ~ 0.4-0.7 keV, plus a power-law (PL) tail extending up to 



10 keV (IMereghetti 



(|2008r ) and references therein). In recent years, a new p ower-law componen t 



extending up to ~ 150 keV was discovered w i th the INTEGRA L satellite ([Kuiper et al 



for recent reviews 



Woods & ThompsonI (l2006h . 



Mereghetti 



(l2006h : 



20081) . This component has a very flat 



uF^ spectrum and is thus clearly distinct from the lower energy PL. A seizable fraction of the 
bolometric emission from SGRs/AXPs is emitted in the 10 - 150 keV energy range, and it is likely 
comparable to that at lower energies, below 10 keV (lMereghettill2008[) . The persistent emission of 
these sources gets temporarily enhanced at bursts (SGRs/AXPs) and Flares (SGRs only), but such 
variations are usually moderate - normally at most by a factor of a few. 

Transient sources, on the other hand, are characterized by a much weaker persistent X-ray 
emission, which in some cases still remains undetected. When detected, it is generally found to 
be at a lower level than rotational energy losses (Lx < |-Erot|)- Transient sources are distinctively 
characterized by X-ray outbursts during which, in addition to the emission of sub-second long 
bursts (similar to the bursts of persistent sources), their persistent X-ray emission is subject to very 
large enhancements, by ~ 2-3 orders of magnitude, making them temporarily as bright as the per- 
sistent sources. In fact, most of these objects have been discovered thanks to such events. In out- 
burst, their luminosities largely exceed rotational energy losses (Lx,outburst ^ l-E mt l ), and gradu 



ally re turn to the much lower quiescent level over timescales of ~ a few years (cfr. iRea & Esposito 



(|201ll) for a recent review). Data for persistent SGRs/AXPs and transient sources were collected 
from the McGill catalogs . 



^ URL: |http://www.physics.mcgill.ca/~pulsar/magnetar/main.htnil 
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Name 


P 
(s) 


P 

(10-11 s/s) 


(10^^ org/s) 


kT 
(kcV) 


(kyr) 


Erotl 
(10'" org/s) 


-Bdip 
(lOl* G) 


a = 


k = 
i-x/IBrot 


/a = 




SGR 1086-20 


7.6022(7) 


75(4) 


1.6 




0.16 


67 


24 


0.00083 


2.4 


0.096 


SGR 0526-66 


8.0544(2) 


3.8(1) 


14 




3.4 


2.9 


5.6 


0.28 


49 


0.32 


SGR 1900+14 


5.19987(7) 


9.2(4) 


0.83-1.3 


0.47(2) 


0.90 


26 


7.0 


0.037 


4.1 


~0.17 


SGR 1627-41 


2.594578(9) 


1.9(4) 


~0.025 




2.2 


43 


2.2 


0.020 


0.058 




SGR 0418+5729 


9.07838827(4) 


<0.0006 


0.00062 


0.67(11) 


> 24000 


<0.00032 


<0.075 


>5050 


>196 


0.000024 


SGR 1833-0832 


7.565408(4) 


0.439(43) 






27 


0.40 


1.8 








IE 1547-5408 


2.06983302(4) 


2.318(5) 


~0.0058 


0.43(4) 


1.4 


100 


2.2 


0.0032 


0.0056 


0.0013 


XTEJ1810-197 


5.5403537(2) 


0.777(3) 


~0.0019 


0.14-0.30 


11 


1.8 


2.1 


0.0092 


0.11 


0.0063 


IE 1048-5937 


6.45207658(54) 


~2.70 


0.054 


0.623(6) 


3.8 


3.9 


4.2 


0.022 


1.4 


0.0028 


IE 2259+586 


6.9789484460(39) 


0.048430(8) 


0.18 


0411(4) 


230 


0.056 


0.59 


225 


320 


0.049 


CXOU J010043.1 -72134 


8.020392(6) 


1.88(8) 


~0.78 


0.38(2) 


6.8 


1.4 


3.9 


0.65 


54 


0.29 


4U0i42+61 


8.68832973(8) 


0.1960(2) 


>0.53 


0.395(5) 


70 


0.12 


1.3 


>40 


>449 


>0.17 


CXOJ164710.2-455216 


10.6107(1) 


0.083(2) 


~0.0044 


0.49(1) 


200 


0.027 


0.95 


1.7 


16.3 


0.0006 


IRXS J170849.0-400910 


10.9990355(6) 


1.945(2) 


~1.9 


0.456(9) 


9.0 


0.57 


4.7 


1.47 


329 


0.34 


IE 1841-045 


11.7750542(1) 


4.1551(14) 


~2.2 


0.44(2) 


4.5 


0.99 


7.1 


0.37 


219 


0.45 


PSR J 1622 - 4950 


4.3261(1) 


1.7(1) 


~0.0063 


~0.4 


4.0 


8.5 


2.8 


0.0064 


0.076 


0.0019 


CXOUJ17I405.7-381031 


3.82535(5) 


6.40(14) 


~2.2 


0.38(8) 


0.95 


45 


5.0 


0.16 


4.9 


0.82 




RXJ 1856 


7.05 


0.003 


0.00017 


0.063 


3700 


0.0034 


0.15 


55 


5.0 


0.084 


RXJ 0720 


8.39 


0.00698(2) 


0.00337 


0.090(5) 


1900 


0.0047 


0.25 


202 


72 


0.40 


RXJ 1605 


6.88 




0.0001 Id^ 1 


0.096 












0.010 


RXJ 0806 


11.37 


0.006 


0.00033 


0.096 


3000 


0.0016 


0.26 


27 


20 


0.030 


RXJ 1308 


10.31 


0.011 


0.00005 lti§ J 


0.102 


1500 


0.0040 


0.34 






0.0036ti§ J 


RXJ 2143 


9.44 


0.004 


>0.00069 


0.100 


3700 


0.0019 


0.20 


^126 


^37 


^0.053 


RXJ 0420 


3.45 


0.00028 


0.00032 


0.044 


2000 


0.027 


0.1 


115.2 


1.18 


0.66 



Table 1. Summary of the salient timing and X-ray specti'al properties of SGRs, AXPs and XDINs. Data are collected from references cited in the 
text. 



In addition to these "classical" magnetar candidates, we consider the properties of X-r ay Dim 



Isolated NSs (XDINs; 



Kaplan 



(12008b; 



Mereghettil (|201 II) . Kaplan & van Kerkwijk 



2011). These 



are a class of currently seven, ~ 10^ yr old X-ray sources, which are considered as neat examples 
of i solated, coolin g NSs. They have stable and nearly purely thermal persistent emission (but see. 



e.g. lTurollall2009[) . typically described as a single BB with kT ~ 0.05 -0.1 keV. Their spin periods 
are in the same range as SGRs/ AXPs and their inferred dipole fields are all > 10^^ G. As such, a 
possible link with magnetar candidates has long been suspected. Data collected for all sources are 
summarized in Table [TJ 

Although timing data are homogeneuos and can be easily compared for all classes above, 
significantly different values of P have been measured at different epochs in SGR 1806-20 and 
1900-1-14, leading to different values of the inferred dipole fields. We show in all plots the max- 
imum and minimum value of the inferred -Bdip, for either object, joined by a dotted line. In the 
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following discu ssion we conside r , how e ver, the lowest value as a m ore likely indication of the 



dipole field (cfr. 
torque changes). 



Thompson et al. 



(l2002h : 



Woods & Thompson! (|2006r) for discussion of magnetic 



Data on the X-ray emission from different classes, on the other hand, must be compared with 
caution, hnportant spectral differences exist and, in some cases, marked temporal variability makes 
the sample much less homogeneous. We adopt the X-ray luminosity in the (2-10) keV range, as 
reported in the McGill catalog, as a measure of the thermal emission from persistent sources, 
which we call Lx. Given their typical BB temperatures, bolometric corrections are expected to be 
~ 2 - 4 and will be ignored. Note that the contribution from the PL component in this energy range 
is not negligible either, which at least partially balances for the neglect of bolometric corrections. 
Also note that GFs might provide a non-negligible contribution to L^ot in SGrJ^I, and that the hard 
tails up to > 150 keV also contribute significantly to the latter. 

The quiescent emission of transients is less clearly understood. In most of them it appears 
dominated by a BB component with temperature (~ 0.4 keV) and an emitting area much smaller 
than the NS surface. The resulting luminosities are generally < 10^'^ erg s^\ and the bolometric 
corrections to their (2 — 10) keV emission are thus similar to those of persistent sources. 

In XTE J1810, on the other hand, the BB component which dominates the quiescent emission 
has a m uch lower temperature , kT ~ 0.18 keV and an emitting area consistent with the NS 



surface ( Bernardini et al. 



20091 ). This implies a total luminosity < 10 erg s , but a relatively 



large bolometric correction to the (2 — 10) keV flux. 

As we see later the "weak field magnetar" SGR 0418-1-5729, with its extreme parameters, 
plays a very important role in constraining various models. Particularly important is its total X-ray 
luminosity, Ltnt - The current observed value of ^ 6 x 10'^^ erg s~^ in the (0.5 — 10) keV range 



(IRea et al 



2010h might well be a post flare remnant emission larger than its real quiescent level. 



which is unknown. It is interesting to use the energy budget of its recent outburst to estimate a 
minimal average power needed. The decaying X-ray flux following the outburst was monitored 
from June 2009 to September 2010 (|Esposito et al.l |2010|) . Results of this monitoring allow to 
estimate a total fluence in the 0.5- 10 keV energy range corresponding tqj Ai? 3 x 10^° erg, 
for a distance of 2 kpc to this source. Even if the outburst recurrence time was Tree ~ 100 yr, 
this would still correspond to an average luminosity of (Loutb) ~ 10^^ erg due only to such 
outbursts, not much below the current upper limit. 



^ This contribution miglit even be dominant in SGR 1806-20. 

^ This corresponds to an average luminosity of ~ 10'^'^ erg during the outburst. 
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Does such a recurrence time match what we know about outbursts from transient sources? If 
Tage is the age of SGR 0418+5729 and Tbirth is the average time between the birth of magnetars 
in our Galaxy, one expects a rate of outbursts of -Routb = 10 Tage,Myr(7birth,kyr T^ecimY^ yi'^^ 
in the Galaxy. The age of SGR 0418+5729 could be as small as 0.1 Myr (based on the upper 
limit on its quiescent X-ray luminosity, as we discuss in next sections) but, on the other hand, 
their actual birth rate is likely to be a few times larger than our adopted lower limit of 1 per 
kyr (IGaensler et al.lll999h . The ratio Tagc/Tbirth is thus not expected to vary much, in any case. 
Comparing with observations, we note that outbursts from 4 different sources were detected, in 
one year of operation, by the Fermi satellite. Two of these sources were previously unknown and 



could not have been detected if they 



lad been further than a fe w kpc, and are thus detectable only 



from a small fraction of our Galaxy (|van der Horst et al.ll2010|) . This matches reasonably well our 



estimate of -Routb, thus supporting our estimate of (Loutb)- 

We finally note that uncertain distance determinations can affect the inferred luminosities. 
For most objects, different estimates agree to within factors < 1.5. Although non-negligible, this 
typically translates to uncerta inties of a factor < 2 on Ly. The AXP IE 1048-5937 represents, 
however, a notable exception. lOurant & van Kerkwijkl (|2006r) find a distance of 9.0 ±1.7 kpc for 
this source, based on a detailed study of reddening of red giant stars in the fie ld of this AXR This i s 



2005h . 



significantly larger than the 2.7 ± 0.1 kpc reported in the McGill catalogue (IGaensler et al. 
The persistent X-ray luminosity of IE 1048-5937 would accordingly increase by a factor ~ 11, 
reaching the value Lx ~ 6 x 10'^^ erg s^^ instead of ~ 5 x 10^^ erg s^^ as reported in the McGill 
catalog. We consider both values for this source. 



3 MAGNETIC DIPOLE BRAKING AND FIELD DECAY 

The energy loss rate of a magnetic dipole rotating in vacuum is: 

i^vac = - — §— sm 9b , (1) 

o C 

where 9^ is the angle between the dipole and rotation axes. 

= BeqR^ = -Bpo\R^ , (2) 

is the magnetic moment, where R is the stellar radius while B^q and Bpoi are the surface mag- 
netic field strengths at the dipole equator and pole, respectively. More realistically, however, a 
magnetized rotating neutron star is surrounded by plasma rather than vacuum. Three dimensional 
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(3D) force-free numerical simulation (|Spitkovskyl2006l) have shown that this slightly increases the 
energy loss rate to: 



^pis = — ^(1 + sin 9b) . 

The evolution of the spin period P is governed by Lpis = —E^., 



d 



(p2) = 2PP 



l + sin^ Ob) 



(l + sin^^B)^/ 



(3) 

-m(l = I{2'nfP/P^, or: 



cq 



(4) 



dt ' " ^ IC^ ' /C3 ' ' ' IC^ 

Values of Bcq for isolated, spinning down NSs are usually estimated through this formula, using 
the measured values of the spin period, P, and its first derivative, P. Generally the formula for 
an orthogonal rotator in vacuum is used (Eq. with 9b = 90°). That expression is formally 
recovered from the more realistic one given in Eq. dB]) by takingly / = 2/3. In what follows we will 
maintain the explicit dependence on / for all quantities but always specialize to the case / = 2/3 
when making numerical estimates. Defining a characteristic spindown age, Tc = P/ (2P), Eq. dH) 
can be rearranged to read: 



cq 



Ic^ 



1/2 



p 

In 



(5) 



We shall now allow the dipole magnetic field to evolve in time in the above equations and 
use its value at the equator as our reference field strength from here on, Bcq{t) = Bdip{t). For 
simplicity, however, the angle 9b and thus also the factor / = 1 + sin^ 9b will be taken to be 
constant in time. The exact solution for the spin period becomes: 



P\t) = P'{to) + f- 



Ic^ 



to 



Blp{t')dt' . 



(6) 



The spindown depends in a critical way on the time-dependence of magnetic dipole energy, i^dip c< 
i?Jip, since this determines the behaviour of the integral on the r ight-hand si de. 
Following the notations and parameterization introduced by 



Colpi et al. 



(l200nh . we write: 



(7) 



-Bdip(i) — -Bdip.i 



(8) 



dt " '^^ T^iB^p) 
where we define the field decay timescale rd(i?dip) = (^-^dip) "^- The solution of Eq. (|7]) is: 
(l + at/rd,i)-i/" (a^O) , 
exp(-t/rd,i) (a = 0) , 

where i?dip,i = -Bdip(t = 0) is the initial dipole field strength, and rd,i = 'rd(i?dip,i) = V^-^dip.i 
is the initial field decay time. Note that a = 0, which corresponds to an exponential field decay 
on the timescale Xd^cxp = 1/^, is the only value of a for which remains constant. We do not 



^ This is just a formal choice, made for comparison with published data, since the angular factor in Eq.|3]can never be smaller than unity. 
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consider a < (for which the field vanishes within a finite time, t = —Td^Ja). Substituting Eq. dH) 
into Eq. ^ gives: 

(^[l-{l + at/T,,)(--'y-] (a ^0,2), 
P\t) = + f^^Bl^.,T,A i [1 - exp(-2t/rd,0] (a = 0) , (9) 

I iln(l + 2t/rd,i) (« = 2), 

where Pj = P{t = 0) is the initial spin period, at the birth of the NS. For ^ a < 2 the spindown 
essentially freezes out at late times (or at very late times for a close to 2) and the rotation period 
P appraoches a constant value. 



r2 D6 



OO I ^ J 



(10) 



Jc3 ^'P'' (2 - a) ■ 

The field decay thus proceeds at a nearly constant spin at such late times. Usually Pqo ^ P^ and 
in this case: 

V2 



P 



/c3(2 - a; 



7.65 s 



1/2 _ '-"^ ^ .1/2 p3 r-l/2 p 



^ -^dip,i ; 



(11) 



where hereafter we use the notatiorj^ Qn = 10" x Q in c.g.s units. 

It is convenient to rewrite Eq. (|7]) in terms of the spindown time, Tc = P/(2P), using the 
relation dBdip/drc = Pjipdrc/dt: 
1 



dTr 



-Bdip 



(12) 



For the initial conditions Tc = Tc^, Pdip = 5dip,i and Xd = rd,i at t = 0, the solution to this equation 



is: 



Tr(B. 



cl,-DdipJ 



Bdip 



Tc i + TdJ In 



2 /„ \ a 



^dip.i 
Bdip 



(« ^ 2) , 

(a = 2) . 



(13) 



Substitution of Eq. ([8]) into this solution yields: 



2-a 



{ [l + (2 - ")^] (l + - (l + e) } (" ^ 0' 2) 



Td.i 

2 



1 + 



2Tc,i 
Td.i 



exp 



(a = 0) 



(14) 



1 + ^ 

T"d,i 



rc4 + ^ln(l + ;f- 



(a = 2) . 

For the specific case of a = the expression for Tc{t) in Eq. (fT4l) can be easily inverted to obtain: 



To avoid too many subscript, the initial value of the dipole field is indicated simply as Bi^n, when it is normalized to the n-th power of 10. 



t(re) 



2 I 1 + ^ 

Td,i 
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(15) 



but this is not possible for a general value of a. Nevertheless, at early times we can write: 

Equations © and (flOl) imply that: 

-1 



(16) 



P 



- 1 



P 



(17) 



Since one expects P,; ^ Poo, then the term in Eqs. (|T3l) -(fT5l) involving the initial conditions (i.e. 
that with Tc^i) can be neglected at t ^ t^^i, where Tc{tc;i <^t <^ i/a) ^ t. 
The limits of the expressions above at late times, t ^ r^^i/a, are: 

(«t/rd,i)-^/" (a>0) , 



Bdip{t > rd,i/a) ^ Pdip,i 



exp(-t/rd,i) (a = 0) 



(18) 



and 



P(t > rd,i/a) 



[a 



(a = 0) , 
(0 < a < 2) 
(a = 2) , 
(a > 2) . 



(19) 



P^{a) [l-iexp(-2t/rd,i)] ^ P, 
Poo(«) [1 - i(«t/rd,i)("-2)/a] _ p^(^) 

Poo(a = 0)v/ln(2t/rd,i) 
Poo(4-a) (at/rd,i)("-2)/2" 

For a > 2 the period grows at late times as P oc jt(°-i)/2" T;\iis reproduces the familiar P oc P^^ 
evolution for spin-down at a constant magnetic field (a — )■ oo). For ^ a < 2, P approaches 
Poo at late times. We note, however, that for < 2 — a ^ 1 the asymptotic spin period Poo 
is approached extremely slowly, where P(t)/Poo ^ g < 1 at t ^ (rd,i/Q;)(l — (yi)^"/^^^") — >■ 
(rd,i/2)(l - e.g. P reaches half of its asymptotic value at t ^ (rd,i/a)2"/(2-a) _^ 

rd,i2^/*^^~°). When Poo is approached, then Eq. dS]) (or equivalently, Eq. |[T3l ) implies that: 



B. 



dip 



1.07 X 10-/-''^/if (^) 



-1/2 



G 



(20) 



Note that, although this expression is independent on a, it holds only for a < 2. This simple 
scaling will indeed hold for any scenario in which the spin period freezes, provided the appropriate 
value of Poo is used. In fact, it corresponds to the general relation [51 with the spin period kept 
constant to the asymptotic value. 

From Eqs. (fT3l) and (fT4l) we obtain the late time behaviour of Tc, at Pdip ^ ^dip.i or t r^X- 
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Figure 1. Different cases of field decay witli asymptotic spin period. Left Panel: Dipole magnetic field, Bdip. as a function of spindown age, Tc, 
for 4 selected values of the decay index, a ^ 2. All curves have the same initial i?dip,i = 10^^ G. Smaller values of the asymptotic period Poo 
con'espond to lower asymptotic curves where B^jp oc Poo / Right Panel: The ratio t/rc between real age and spindown age as a function of 
the spindown age, for the same models. The y-axis gives the correction factor that must be applied to the measured spindown age, Tc, to derive the 
real age of an object. Lower values of a con'espond to much younger objects, at a given Tc, because they produce a faster decay of the field and, 
accordingly, of P. 



2 



2-a 



B, 



dip 



-Bdip 



a~2 



^dip.i 
^^dip 



-Bdip.i 
Bdip 



Td.i / at_ 
2~a \ rd,i 



Td.i 



2/a 



a— 2 



(« = 0) , 
(0 < a < 2) 
(« = 2) , 
(« > 2) , 



(21) 



The are summarized in Fig. [T] summarizes the main results of this section. It shows the expected 
relation _Bdip vs. Tc and the evolution of the ratio {t/rc) with Tc, for 4 different values of « ^ 2. 
All plots are obtained for the same initial value -Bdip.i = 10^^ G and assume a normalization of the 
decay time, rd,i = 10^yrs/-Bf]^5 yrs. 



4 SUMMARY OF PHYSICAL MECHANISMS CAUSING FIELD DECAY 



Gene ral modes of fled decay in non-superfluid NS interiors were studied by lGoldreich & Reiseneggei 
19921 (GR92 hereafter), who identified three avenues for field evolution: Ohmic decay, ambipolar 
diffusion and Hall drift. While the first two mechanisms are intrinsically dissipative, leading di- 
rectly to a decrease in field energy, the third one is not. GR92 proved its potent ial relevance , 



however, for the transport and dissipation of magnetic energy within NS crusts (cfr. 
by speeding up ohmic dissipation of the field. 



Jones 



1988h . 
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The analysis by GR92 showed that ambipolar diffusion is conveniently split in two different 
components, according to their effect on the stable stratification of NS core matter. The solenoidal 
component does not perturb chemical equilibrium of particle species, thus its evolution is opposed 
only by particle collisions. The irrotational mode does perturb chemical equilibrium, so it also 
activates /3-reactions in the NS core. At high temperatures (T > 7 x 10*^ K) the two modes are de- 
generate, as /3-reactions are very efficient and particle collisions represent the only effective force 
against which ambipolar diffusion works (GR92, TD96) The distinction between the two compo- 
nents becomes essential at lower temperature, as perturbations of the chemical equilibrium are not 
erased quickly, thus significantly slowing down the irrotational mode (GR92, TD96). Following 
GR92, the relevant timescales for the two modes of ambipolar diffusion can be written as: 



-^t^ ~ 3 X 10^ ^TT^ yr 



(ir) _ 5 X 109 (,) 

'ambi — 2"6^2 J ^ ™,bi • 

8 15 

Note the strikingly different dependence on temperature of the two modes, which turns out to be 
a key factor. Since field dissipation releases heat in the otherwise cooling core, a balance between 
field-decay heating and neutrino cooling (through modified URCA reactions) is expected to be 
reached. This determines an equilibrium relation between core temperature and str ength of the core 



magn etic field. The two following relations, for either mode, are obtained (TD96; 



DairOsso etal. 



2009L hereafter DSS09): 



T^^) ~ 2 7 B^'^ ( ^\ 

^8,eq - ^-'^15 [q^ J 



-2/3 



where pis is the density of matter normalized to 10^^ g cm^^. Using these relations we eventually 
obtain expressions for the relavant timescales as a function of B alone: 

i'bi - ^"g«/5° (gf) (a = 6/5), 



Ci, - 3.7xl0»(^j [^^^J („ = 14). (24) 

Hall-driven evolution of the magnetic field is characterized by the timescale THaii = A'KeniL?' / B . 
Here Ue is the electron number density and L is a characteristic length on which significant gradi- 
ents Tie and B develop. GR92 argued that the main effect of the Hall term would be that of driving 
a cascade of magnetic energy from the large scale structure of the field to increasingly smalle r 



scales. The ohmic dissipation timescale is Tohm = ^t^OqL^ jc? oc (GR92. ICumming et al. 



20041) . 
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where (Tq is the electrical conductivity of NS matter. Thus, very efficient dissipation of sufficiently 
small-scale structures would eventually be reached. Such a "turbulent" field evolution would be 
of particular relevance in NS crusts, where (1) matter density is lower than in the core and the 
Hall timescale is then short enough (2) ions are locked in the crystalline lattice and B-field evolu- 
tion is thus coupled to the electron flow only. The importance of the Hall term relative to ohmic 
dissipation is typically quantified by the Hall parameter, uo-qt = eBT/{m^c) = Tohm/THaii- Here 
m^: = Ep / ^ rrie is the electron effective mass, E-p being its Fermi energy, and r is a typical 



electron collision time. iCumming et al.l (|2004) have shown that the Hall parameter in a NS crust is 
always much larger than unity, if the magnetic field is > 10^^ G, apart from, possibly, the lowest 
density regions of the outer crust (cfr. their Fig. 4). For weaker fields, on the other hand, the ohmic 
term might largely dominate at temperatures lower than a few x 10^ K. Hence, a prominent role of 
the Hall term is expected f or magnetar- strength fields. 



Gumming et al. 



([200^ highlighted the prominent role played by a realistic electron density 
profile in determining the Hall evolution of the crustal field. As Hall modes first get excited at 
the base of the crust, where p ^ 10^^ g cm^^ and the Hall time is longer, they can propagate to 
the surface with their wave-vector k progressively decreasing, because of the decreasing electron 
density. The ohmic dissipation time decreases accordingly and the overall decay rate of the field 
is set by the longest times cale at the base of crus t. Using the local pressure scale height P/ (pg) as 
the natural lengthscale L. iGumming et al.l (|2004j) derive the expression: 



THall 



7/3 

~ 1.2 X yrs 

-Bdip,15 



(25) 



for the Hall decay timescale in NS crusts. Note the different normalization of the density compared 
to ambipolar diffusi on timescales, reflecting the different locations of the t wo processes. 



Several authors (Vainshtein et alj 



2000 



Rheinhardt & Geppert 



20021: 



Pons & Geppert 



201C) 



further investigated this scenario numerically, generally confirming the tendency of crustal fields 
to develop shorter scale structures on the Hall timescale. However, numerical calculations may not 
fully support the idea that the Hall timesca l e actually characterizes the decay timescale of mag 



netic modes (see 



Kojima & Kisaka 



Shalvbkov & Urpin 



1997 



HoUerbach & RUdiger 



2002 



Pons & Geppert 



20071 : 



20121). The situation in this case is likely more complex than the basic picture 



given above. 

Alternatively, ohmic dissipation can proceed at a fast rate if the electric currents and the field 
are initially roo ted in relatively outer lay ers of the crust, where the electrical conductivity cTq is 



rather low (cfr. 



Pethick& Sahrlinglll995h . Due to the subsequent diffusion of electrical currents 
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into deeper crustal layers, (Tq progressively increases and field decay slows down accordingly. The 
decay of the dipole field can effectively be described as a seque nce of exponentials with a growing 



characteristic time rH,exp(t). In this framework it can be sho\y n (|Urpin. Chanmugam & Sang 



Urpin. Konenkov & Urpin 



19971 : 



19941 : 



Urpin & Konenkovll 2008h that, although the ohmic dissipation 



rate is field-independent, the resulting field decay proceeds as a power-law in time, after a short 
"plateau" which is set by the initial distribution of currents. The decay index of the power-law is 
determined by the rate at which cTq changes with time. This depend s on both the depth reached 
by cu rrents, h{t), and the temperature of the crust, Tc(t). It is found (lUrpin. Chanmugam & Sang 



1994) that a power-law decay with index =3/2 results when (Tq is independent on temperature, 
as is the case when impurity scattering dominates the conductivity. An index = 11/6 ~ 1.83 
is obtained, instead, wh en is dominated by elec t ron-p honon scattering and, thus, scales with 



T-^ (cfr. Fig. 2 and 3 of 



Gumming et al 



Urpin. Konenkov & Urpin ( 1997|) and the discussion of ohmic decay in 



20041) . 

Note that an asymptotic spin period Pqo ^JT^i 5dip,i exists also in this model. Since the 
decay time rd,i is field independent, the expected scaling between Poo and Bdip,i is linear, like in 
the exponential case (a = 0) previously discussed. 



5 OBSERVATIONAL EVIDENCE FOR FIELD DECAY IN NEUTRON STARS WITH 
STRONG MAGNETIC FIELDS 

We begin by assessing whether the distribution of dipole magnetic fields of magnetar candidates, 
as inferred from timing observations, can provide indications for field decay. Fig. [2] shows the 
inferred dipole fields, Pdip^ of magnetar candidates plotted versus their spindown age, r^. This 
represents an alternative projection of the usual P — P diagram. 

Two points stem most clearly from Fig. |2l First, the absence of old objects (with relatively 
large spin-down ages) with strong dipole fields and the tendency for objects with increasingly 
stronger fields to be found only at increasingly younger spin-down ages. The dashed line in Fig. |2] 
represents the scaling of Eq. (l20l) . expected from field decay with a < 2, using Poo = 11-77 s that 
corresponds to the longest measured spin period (for AXP 1841-045). The apparently prohibited 
region in parameter space is equivalent to the existence of a limiting spin period for the objects 
considered. As discussed in previous sections, precisely such an asymptotic spin period is expected 
if the dipole field decays and its decay is governed by a physical mechanism with a <2. 

Second, all "old" sources (spin-down age Tc > 10 kyr) lie in a narrow strip corresponding to 
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Figure 2. Inferred Dipole Field (B^ip = 3.2 X 10^^ V PP G) vs. characteristic spindown age (tc = (1/2)P/P) for SGRs, AXPs, transients and 
XDINs. The dashed line in Fig.[2]represents the scaling of Eq. )20K expected from field decay with q < 2, using Poo = 11.77 s that corresponds 
to the longest measured spin period (for AXP 1841-045). The magenta arrow indicates the current upper Umit on the position of SGR 0418+5829. 
The red triangles joined by dotted lines represent the whole range of values spanned by SGR 1806-20 and SGR 1900-1-14, as explained in the text 
(cfr. Tab.[T). 

3.5 s < Poo ^ 11 s, i.e. within a factor of ^ 3 in their spin periocj^. In the context of magnetic field 
decay models, the fact that this strip is so narrow can be interpreted in one of two ways. Either 
(1) all "old" AXPs, XDINs and SGR0418+5729 came from a narrow distribution of initial dipole 
fields -Bdip.i, in which case a relatively large range of a < 2 values would be compatible with 
observations, or (2) a is sufficiently close to 2 that field decay has largely washed out the spread 
in -Bdip,i values, which could have been significantly larger in this case. 

The distribution of sources shown in Fig. [21 can be understood as follows. For young objects, 
t <^ Td^, the magnetic field does not have time to decay and it is almost constant. This implies a 
constant B and a constant Eb^^p oc -Bdip"- ^^^^ stage the period P and the spin-down age grow 
while the dipole field B^ip doesn't vary (corresponding to a horizontal trajectory in the i?dip-Tc 
plane). Moreover, the object has a nearly constant power output due to magnetic field decay up to 
Tc < Td i, at which point the power drops following the decay of the magnetic field. This evolution 
scenario implies that objects with age t ^ rd,i are most likely to be detected close to rd,i since, for 
a constant power output, they spend more time at Tc ~ Td^. 

This is particularly expected for SGRs, since they are detected predominantly through their 
bursting activity, which is thought to be directly powered by the decay of their magnetic fielcj^, 

^ Note that, with the notable exception of the XDIN RX J0420, all other sources would lie between 7s and 11 s,arange~ 1.6 in spin period. 
* This assumes that their bursting activity is powered by their dipole field, rather than by their internal field. If the latter decays on a longer 

timescale this might account for SGRs further along the dashed line in Fig. [2] 



From SGRs to "weak field magnetars " 17 



« = (Td,exp = yrs) a = 1, Tj = (lO^/Bjip ,5) yrs 
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Figure 3. Different models of field decay with asymptotic spin period, (a < 2) compared with the inferred dipole values of magnetar candidates 
from timing measurements. The magenta arrow represents the current upper limit to the position of SGR 0418+5729. The plotted curves don't 
represent fits to data but the normalization of the decay timescale (A) was chosen, for each value of a, to match the position of sources. 



rather than through their quiescent emission (which might have a non-negligible contribution from 
the NS residual heat). Therefore, the fact that we do not detect any SGRs with magnetic fields of 
several x 10^'^ G but larger spin periods (P ^ 10 s), strongly supports a decay of the dipole field 
on a timescale of rd,i ^ 10^(i?dip,i/10^^ G)"^ yr, in these objects. Note that, if the SGR's magnetic 
dipole field did not decay on such a timescale, a large population of SGRs with spin periods of tens 
or even hundreds of seconds would exist. Such large period SGRs should be easy to detect if they 
maintained similar magnetic fields and thus similar bursting activity as the observed SGR sample. 
Even objects with much weaker or no bursting activity but larger spin-down ages are detected (e.g., 
AXPs, transients, XDINs) which, if there was no magnetic field decay, would imply similarly 
larger ages. Therefore, we find it highly unlikely that long period SGRs, with P ^ 10 s, exist 
in much larger numbers than the observed SGR sample but are not detected for some unknown 
reason. 

Objects with ages t <^ Td^, which populate the lower-left region of the B^-Tc plane, are 
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Figure 4. Different models of field decay with asymptotic spin period {a < 2) in the i?dip vs. P plane. The magenta arrow indicates the current 
upper limit to the position of SGR 0418+5729. The curves correspond to the same models of the previous figure. 



detected with proportionally smaller probability, unless their bursting/flaring activity is stronger at 
younger ages, resulting in a greater detectability that could compensate for the short time spend 
in this part of the diagram. Relatively old objects with an age t ^ rd,i will have reached their 

— 1/2 

asymptotic period Poo (for a < 2) and will thus be found on the asymptotic line, Tc . 

These will have a much lower luminosity than their younger brethren and will be detected only if, 
e.g., they are sufficiently close to us or if they maintain some level of bursting activity. 

Finally, we note that 5 out of 7 transient SGRs/ AXPs lie, in the Sjip vs. t^. diagram, below the 
asymptotic line. This suggests that Tc is not much longer than their rd,i. CXO J164710. 2-455216 
and, of course, of SGR 0418+5729 represent two notable exceptions, as both appear to lie well 
on the asymptotic line (hence, Tc > rd,i). In particular, the former object has both Tc and Pdip 
comparable to the persistent AXP 1E2259+586. 

We summarize the effect of field decay in Figs. [3] and IH We also show in Fig. [5] the corre- 
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Figure 5. Tracks in the P vs. P plane expected for the different cases of field decay previously discussed, compared to measured source positions. 
Plotted curves do not represent fits to data. Different symbols and lines are as in previous figures. 



spending tracks in the usual P-P diagram. As already stated, previous plots just represent different 
projections of these fundamental measured quantities. 

We can estimate a minimal spread in i?dip,i from Fig.[2]by restricting attention to objects lying 
well below the asymptotic line, so that t^. < r^^. The strongest value of the field is 8 x 10^"^ G 
in SGR 1806-20, while the weakest is 1.8 x 10^^ G in SGR 1833-0832. These numbers imply a 
minimal spread ~ 4.5 for the initial field distribution of the whole population. 

For a given value of a < 2 (and a fixed normalization A), Eq. (fTTj) can be used to derive a 
general relation between the observed spread in Pqo at late times and the initial spread in -Bdip,i 
that produced it. For a distribution of initial dipoles in the range -Bi,min < -Sdip,i < -Bi,max (or 

1 ^ ^dip,i/^i,min ^ — -Si,max/-^i,min) ^ Spread .Poo,min ^ ^oo ^ -Poo,max (O^ 1 ^ Pco/ Poo,mm ^ 

Ap = Poo,max/-Poo,min) is cxpcctcd. Therefore 

= A?/(^-") ^ log Ap ^ log = ^ log (^) ^ ^ log Ap .(26) 

V ^oo,min / ^ V -^i,min / 

In order to quantify Ap from the data, we evaluate the average spin period, (P), and standard 
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deviation, crp, of those sources which are old enough to be considered close to the asymptotic 
line, e.g. having Tc ^ yrs. The resulting subsample contains 13 sources, with (P) ^ 8.42 
s and CTp ~ 2.2 s. The parent population would thus be characterized by a 2a spread of ~ 3.3, 
corresponding to the ratio between the slowest and fastest spinning objects of the subsample. We 
plot relation (|26l) in the left panel of Fig. [6l where the two curves correspond to the la and 2a 
values of Ap, as derived above. 

A totally independent constraint on a can be derived considering the very weak, thermal, X - 



Rea et al 



(12010) ') 



ray emission of SGR 0418+5729. In relative quiescence (on 2010 July 23rd, 
this is dominated by a black-body component with a temperature of kT = 0.67 ±0.11 keV, 
and corresponding 0.5 - 10 keV luminosity, Lx ~ 6.2 x 10'^^ erg s"^ In the earlier, more active 
period (June to November 2009), the temperature gradually decrease d from kT ~ 1.0 ke V to ~ 



0.8 keV and Lx was between one and two orders of magnitude higher (lEs posito et al.l 



2010D, while 



the corresponding emitting areas gradually decreased. This suggests that even this weak X-ray 
emission is more likely to originate from a localized heating event associated to the recent bursting 
activity, rather than being powered by the secular cooling of the NS. Hence, the luminosity of 6 x 
Iq3i gj,g g-1 considered as a solid upper limit to the quiescent X-ray emission of this object. 

This value c ompares well with the X-ra y luminosity of ~ 10^ yr-old, passively cooling objects, like 



B0656-I-14 ( Yakovlev&Pethick 



20041 ). while younger 10'^ - 10'' yr old isolated, passively cooling 



NSs are at least one order of magnitude brighter than the upper limit for SGR 4018+5729. Since 
field decay is likely to provide additional heat in the latter object, the age of B0656+14 represent 
as a robust lower limit to the true age of SGR 0418+5729. 

In the right panel of Fig.[6]we plot the estimated age of SGR 0418+5729 as a function of a, for 
the whole range of values that give an asymptotic spin period (0 ^ a < 2). The line corresponding 
to the above lower limit, 10^ yr, is drawn for clarity. Note that as a approaches 2 from below the 
true age of SGR 0418+5729 becomes progressively closer to its spindown age. 

Combining the constraints from these two figures, we can rule out, or at least to consider very 
unlikely, values of a < 1. Combined with the earlier constrains this implies 1 < a < 2. 

The points discussed above are illustrated quantitatively in Fig. |3l where trajectories in the 
^dip - Tc plane are plotted for different values of a and for a given initial distribution of dipole 
fields, i?dip,i, comparing them to the objects in our sample. Fig. |4] shows trajectories in the -Bdip- 
P plane for the same sources, same field decay models and i?dip,i distributions as Fig. [3l This 
represents an alternative way of illustrating the same argument. Below we draw some conclusions 
from these plots regarding the viability of different a values. 
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Exponential decay (a = 0) can explain the observed distribution of sources in the i?dip - Tc 
or i?dip - P planes if the field decay timescale is particularly short, = t^^i ^ 1 kyr. Longer 
timescales would fail to halt the spindown at sufficiently short spin period, so this is quite a strict 
requirement. This model could work if all the sources considered here came from a very narrow 
distribution of initial dipole fields (as can be seen in Figs. |3] and Hj). In particular, Fig.|4]shows that 
even a factor 4 spread in -Bdip,i would lead to a significantly wider distribution of spin periods than 
observed (i.e a factor of 4 spread in Pqo compared to the observed factor of < 3). Indeed, Eq. (fTTI) 
shows that the value of Pqo is linear in the initial field, B^ ^, for a = 0. 

Given the short timescale required by the exponential decay, the implied age of SGR041 8+5729 
would have to be younger than, at most, 5 kyr, as can be derived from Eq. (1151) . This is in sharp 
contrast with the minimal age that we derived for SGR 0418+5729 based on its weak thermal 
emission. Additionally, the age of SGR 0418+5729 would not be much larger than that of other 
transients with much smaller t^, while its X-ray emission would be a factor 10-50 lower. Ex- 
plaining such a fast decrease in Lx at these young ages would also represent a major challenge. 
Overall, exponential decay of the field appears ruled out based on the very weak emission of SGR 
0418+5729. 

The case a = 1 provides a reasonable description of the distribution of sources in Fig. [3] and 
Fig. m Although this case suggests a straightforward relation to the basic Hall decay mode (cfr. 
sec. HI), note that the curves in Fig. [3] and |4] were drawn by adjus ting the normalization of the decay 



timescale. A, to be 10 times smaller than the value provided by 
Namely, we assumed 



Gumming et al. 



(l2004h and GR92. 



10^ 

Td ~ yrs (27) 

dip, 15 

The decay for sources whose initial field was larger than ^ 2 x 10^'' G would be to slow with the 
Hall decay time of Eq.|25]and in this case most AXPs/SGRs would evolve to significantly longer 
spin periods at later times. These older counterparts would occupy a region right above XDINS 
and SGR0418+5729 in the (tc, B) plane (Fig.O where no object is actually found. Note, however, 
that the timescale of Eq. [25] could match our empirical scaling if the decay of Hall modes were 
regulated by processes occurring at somewhat lower crustal densities, < 3 x 10^'^ g cm^^, well 
above the crust/core interface. 

This model for field decay implies that all sources come from a distribution of initial dipoles 
in the range (0.2 — 2) x 10^^ G, although most of them (20 out of 23) were between 5 x 10^^ 
G and 2 x 10^^ G. The XDIN RX J0420 would represent a notable exception, having reached a 
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Allowed spread in Bjjp | vs. a 



Ageof SGR0418+5729 (0<a<2) 



Ap = 1 .75 
Ap= 3.3 ■ 
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Figure 6. Constraints on the value of a. Left Panel : The allowed spread in the initial magnetic field distribution as a function of a, for 3 different 
values of the spread in asymtpotic spin period. Right Panel : The inferred age of SRG 0418+5729, as a function of the chosen value of a. The likely 
lower limit ~ 10^ yrs to the age of this source is drawn as a horizontal line. 



remarkably short asymptotic spin period of ~ 3.45 s because its initial dipole was weaker. The 
two transients IE 1547-5408 and SGR 1627-41, which are slightly below the range of the other 20 
objects, would also populate the weak- field tail of the distribution of -Bdip.i- 

Overall, this scenario appears to account well for the observed distribution of sources, although 
a full self-consistent population synthesis model is needed to verify this quantitatively. 

Phenomenological decay law with intermediate value of 1 < a < 2 (e.g. a ^ 1.5). Although 
there are no existing models for field decay predicting 1 < a < 2, we choose this particular value 
of a as representative of cases in which a wider spread in i?dip,i (^b ^ 1) is allowed, despite the 
observed modest spread in Poo (Ap ^ 3). 

Fig.[3]shows i?dip vs. Tc trajectories with this choice of a. One can see that these converge to a 
narrow strip despite an initial wide distribution of -Bdip,i- This scenario can well reproduce also the 
distribution of sources in the i?dip - P plane, including the apparent absence of spin periods longer 
than c:^ 12 s. Note that the normalization for the decay timescale, rd,i = 10'^ yr for i?dip,i = 10^^ G, 
was chosen to match the scaling of Eq.[20l represented by the dashed line in Fig. [21 

An effective power-law decay w ith index ^ 1.5 — 1.8 is expected in the ohmic decay model by 
Urpin. Chanmugam & SangI (|1994|) . This is not completely equivalent to our phenomenological 
decay law, though, since in that model the ratio -Bdip(i)/-Bdip,i is still a universal function, deter- 
mined by the field-independent parameter rd,i. In particular, this implies a linear relation between 
Poo and Pdip,i- Thus, in order for the asymptotic spins of our sources to all fall in the observed 
narrow range, a correspondingly narrow range in Bdip,i is required. This is the same problem as 
in the a = case. To circumvent it, one has to assume that Td^, the initial decay time, is itself a 
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function of i?dip,i- This extra assumption would make the crustal decay model totally equivalent to 
our phenomenological model. Since rd,i is determined by the initial location of the electric currents 
that sustain the field, a strict (anti)correlation between the initial field strength and the initial depth 
at which currents flow is implied. This is far from trivial and an account of its implications is left 
for future study. 

The limiting case, a = 2 is shown for the sake of completeness. As the previous case, it pro- 
vides a narrow range of asymptotic spin periods, even starting from a very wide distribution of 
initial magnetic field values. However, as already discussed, it does not have an actual asymptotic 
spin period, which implies that spin periods of order 20 s would be expected in older objects. 
Thus, this case is disfavored by the data compared to a ~ 1.5. Also in this case, the normalization 
for the decay timescale was chosen arbitrarily as above. 

It is also possible that two different mechanisms for field decay operate, respectively above 
or below some threshold value of the dipole field, B^. We refer to these as the early and late 
mechanism for field decay, according to the time when they dominate, and denote their decay 
index as Ue or a, respectively. At each value of i?dip the mechanism with the shorter decay time 
determines the overall field decay and, of course, the two times are equal at B^, with the value r^,. 

In such a scenario, a wide distribution of -Bdip,i at birth can also result in a narrow distribution 
of Pr^. As such, it is only meaningful for a much below 2 (say a = \) and significantly larger 
tharo 2. Sources with i?dip > 5* would initially evolve along a shallow trajectory, oc (t/Td,i)~^/"'' 
in the i?dip - Tc plane. Because > 2, trajectories for different initial fields would converge to a 
narrow bundle of curves while Tc < r*. Upon reaching the point where = r^, they would then 
meet the line corresponding to the late mechanism (a < 2) and follow that beyond . 

Applying this reasoning to Fig. [2] we derive the values of 5* ~ 6.9 x 10^'' G and r* ~ 4.5 kyr, 
corresponding to the position of the AXP IE 1841-045. However, only SGR 1806-20 is above this 
threshold, implying that the early mechanism could still be dominant only in it. Stated differently, 
in all the observed sources - but one at most - invoking a second, early decay mechanism doesn't 
help to reduce the initial spread in i?dip,i- 



If Oe < 2. it will cause sources to reach the asymptotic spin period before the late mechanism becomes operative, thus correspondingly to only 
one effective mechanism. 
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6 THE PERSISTENT X-RAY LUMINOSITY 

Having assessed the decay of the dipole field in magnetar candidates, we now turn to test whether 
such a decay can account for their observed X-ray luminosity, Lx, which typically exceeds their 
spin-down power, | E^ox. \ ■ 

We define the magnetic luminosity of the dipole field, Le.dip = 2i?dip/''"d» as the available 
power of a dipole field decaying on the timescale r^. Here, the total energy of the field is Edip = 
{AtcRI/3)BI.J{8tt) = RlBljQ ~ 1.7 x lO^^E^jp^gi?^^ erg and, following the definition of 
Eq.©: 

- _ _:^^(^dip) _ RlBlp _ 2^dip _ . o ^ in^^ r2 dS -i -i .2o^ 

^B,d,p- - 6 - - -5.dxiU i^dip,i5^*,6^d,kyr erg s .(28) 

Fig. |7] depicts a comparision of the total available magnetic power for the same four field 
decay models examined in sec. [5] with the observed X-ray emission of SGRs, AXPs, transients 
and XDINs. The three main groups of sources (SGRs/AXPs, transients and XDINs) populate 
three separate regions in parameter space and in the following we comment on them separately. 
We focus, first, on the persistent AXPs/SGRs, for which results and conclusions are much clearer, 
and comment later on the more problematic XDINs and transients. 

6.1 Persistent SGRs, AXPs and SGR 0418+5729 

Overall, the evolution of Le.dip with as derived from the above models for field decay does 
not match well the distribution of sources in parameter space. In particular, relatively old sources 
(Tc ^ 10^ yrs) tend to be more luminous, in the 2-10 keV energy range, than the maximum 
available power from decay of the dipole field namely, Lx > -^^B.dip- These old sources thus 
provide the strongest indication that even the decay of their dipole field is not able to explain the 
X-ray luminosity. 

An exponential decay with a = and with a short appears consistent with all observations. 
However, as previously discussed, an exponential field decay on a ~ 1 kyr timescale would imply 
an implausibly young age 5 x 10^ yrs) for SGR 0418-1-5729 and would require a very narrow 
initial magnetic field distribution. On these grounds, exponential decay of the field is discarded 
anyway. 

Decay with a = 1 falls short of the observed emission of the two AXPs with Tc > 10^ yr. 
Although the mismatch is apparently by a small factor, we stress that curves represent the total 
available magnetic power, while data points include only the observed emission in the 2-10 keV 
energy range. The bolometric luminosity is a few times larger. Gravitational redshift also intro- 
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Figure 7. A comparison of tlie maximal bolometric luminosity as a function of Tc expected from the decay of the dipole field, according to the 
models of section 5, with the observed X-ray luminosity of different classes of magnetars. The magenta arrows indicate current upper limits to the 
spindown age and Lx of SGR 0418+5729. 



duces a non-negligible correction to the observed luminosity (see sec. 17.11) . Finally, the radiative 
efficiency may well be less than unity since field decay is also expected to dissipate energy in other 
channels (such as neutrino emission or bursting activity). The presence of two AXPs beyond the 
upper curve must therefore be considered as a significant failure of this model. 

For a = 1, the observed X- ray luminosity of SGR 0418+5729 is 30 times larger than the 
available power from dipole field decay, according to Eq. [28l Although, as discussed in section 
[J2l this may overestimate the quiescent power, we don't expect the average power output to be 
lower by more than a factor of 10 than that and, even with this correction, the dipole energy is still 
insufficient. 

For values of a > 1 the mismatch between curves and observations is striking. The distribution 
of persistent sources in the Lx - Tc plane is at best marginally consistent with being powered by the 
decay of the dipole field, if a = 1, while it is totally inconsistent for larger values of a. In general, 
the X-ray luminosity appears to decay on a longer timescale, and to be associated with a larger 
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energy reservoir, than can be provided by the decaying dipole fields. Sources are indeed found at 
a fairly constant level of X-ray luminosity up to Tc ~ 10 kyr. The few older objects known clearly 
show a decline of Lx with (spindown) age which, however, is flatter than expected from the decay 
of the dipole field. 

In the framework of the magnetar model, the most natural explanation for these findings is 
that the persistent X-ray emission of these sources is powered by the decay of an even larger 
field component, stored in their interior In the next section we try to assess the location of such a 
(presumably magnetic) energy reservoir and its most likely decay modes. 



6.2 X-ray Dim Isolated Neutron Stars 

The thermal X-ray luminosity of XDINs is much dimmer than that of persistent magnetar candi- 
dates. It is even possible, in principle, to explain it within the so-called "minimal cooling" sce- 
nario for passively cooling NSs. However, it is difficult to reconcile the relatively bright emis- 
sion of RXJ 0720 with its apparent old age (~ 2 x 10 ^ yr) within that scenario. This suggests 



a possible role of strong internal heating in this object (|Page et al. 



20041) . On the other hand, the 



possibility of mild internal heating in most XDINs has also been recently suggested, based on a 
tend ency for their effective temper atures to be higher than that of normal pulsars with the same 



age ([Kaplan & van Kerkwijkll201 Ih . 

From Fig. |2] we conclude that XDINs should be interpreted as objects whose dipole field has 
decayed substantially, hence their true ages must be younger than their spin-down ages, which are 
narrowly clustered around Tc ~ a few x 10^ yr. This favors passive cooling models in accounting 
for the effective temperatures and X-ray luminosities of XDINs, possibly removing the need for 
internal heating in all these sources. On the other hand, for a = 1 the energy released by the dipole 
field decay is not negligible compared with their X-ray luminosities (see Fig.|7]). Whether this can 
significantly affect their temperatures as compared to normal pulsars can only be checked through 
detailed cooling modeling, which are beyond the scope of this work. 

If, however, detailed modeling will reveal that these luminosities are too large for passive 
cooling then we have to consider heating sources. RXJ 0720 stands here as a unique object. Even 
with a = 1, Eq.[28]yields LB,dip ~ 8 x 10^^ erg , which is a factor ~ 5 lower than its measured 
Lx- This could not be explained by L^Aip and would thus hint to the presence of an additional 
energy reservoir, siimilar to magnetar candidates. In this case, some level of bursting activity would 
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Name 


Tc [kyr] 


tagc [kyr] 
(a = 1) 


tag 

(« 


;e [kyr] 
= 1.5) 


fcTeft [eV] 


RXJ 1308 


1500 


30 




100 


102 


RXJ 0806 


3000 


40 




160 


96 


RXJ 0720 


1900 


42 




160 


90 


RXJ 2143 


3700 


50 




200 


100 


RXJ 1856 


3700 


66 




350 


63 


RXJ 0420 


2000 


100 




580 


44 



Table 2. Spindown age, true age (for two dilferent values of the decay index, a) and effetive temperature of the six XDINs considered in this work. 



accordingly be expected from this source. For a > 1, the available L^^Aip is insufficient even to 
account for the X-ray luminosity of RXJ 2143 and it is marginal compared with other XDINs. 

We conclude that the decay of their dipole field implies that XDINs are younger than their 
spindown age, Tc, which could improve the match between cooling models and their X-ray emis- 
sion properties. On the other hand, there is no compelling evidence for a dominant contribution to 
their X-ray emission from the decay of the dipole field. The case a = 1 is, at most, marginally 
consistent with this hypothesis while, as a becomes > 1, it is increasingly hard for LB,dip to match 
the observed Lx. 

We can use the observed emission properties of XDINs to further constrain the possible values 
of a. It was recently shown that the spindown ages of XDINs are s ystematically longer, by a facto r 



10, than the ages of NSs with similar effective temperatures ([Kaplan & van Kerkwijk 



2Q11). 



The simplest intepretation of this fact is that, due to dipole field deacy, their overestimate their 
true ages by approximately one order of magnitude. We can compare quantitatively this statement 
with models for decay of the dipole field. For a given value of a, the formulae of §|5]allow to derive 
the value of -Bdip,i corresponding to each individual XDIN from its and i^dip (or equivalently, P 
and P) and, from its Tc, obtain its real age, tagc- These estimates, for the cases a = 1 and a = 1.5, 
are summarized in Table [2l For the sake of completeness, we report in the last column of this 
table the effective temperature, kT^s, for each source, as derived by the spectral fits. Note that the 
derived age distributions match fairly well the distribution of effective temperatures, as opposed 
to the Tc distribution. It is clear, however, that low a values give corrections to the spindown ages 
of XDINs by a factor < a hundred, which it too large. The case a = 1.5, on the other hand, gives 
just the right correction factors of order ten. 

Finally we note the striking clustering of XDINs at Tc E [1.5 — 3.7] Myr. Together with their 
comparatively wider spread in X-ray luminosities, [< lO"'^ — 3 x lO'^^] erg s~\ it suggests that 
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these sources may have reached some threshold age, at which cooling becomes very efficient, the 
luminosity drops sharply and the sources become unde tectable. This can happen, for example ,if 
the sources enter the photon-cooling dominated regime (|Page et al.1 (|2004j) and references therein). 
In this framework, the actual age of XDINs should be just slightly larger than the time at which 
the transition to photon-cooling occurs. 

From Tab.[2lwe see that, for a = 1, all XDINs would have ages ^ 10^ yr, most of them being 
significantly younger than that. In this case, a transition to photon-cooling at tagc < 3 x 10^ yrs 
would ha ve to be invoked, however, in most "standar d" cooling models this transition occurs at 



> 10^ yr dYakovlev & Peth 



transition at ~ 3 x 10'^ yr (IPage et al 



ck 



2004: 



Page et al. 



2004, 



2004|) . Only the fastest cooling models have this 



201 1|) and are thus just marginally consistent. Even 



in this case, however, the weak X-ray emission of XDINs at such young ages would represent 
by itself a major problem. No cooling model predicts X-ray luminosities < 10^^ erg s^^ at ages 
tage < 10^ yrs. As a conclusion, the case « = 1 is clearly very problematic from this point of view. 

For a = 1.5, on the other hand, all ages are ~ [1 — a few] x 10^ yrs. This allows an easier 
interpretation of the XDINs as cooling NSs having just passed the transition to photon-cooling, 
the older ones being progressively cooler and dimmer. 

The observed properties of XDINs thus strongly argue against a value of a = 1 and are much 
more consistent with 1.5 < a < 2. In particular, the value a ^ 1.5 provides an overall good 
agreement with the main observed properties of these sources. Altogether, taking into account the 
requirement of a sufficiently below 2 to account for lack of periods well above 10 s, we conclude 
that 1.5<a;<1.8is strongly favored by the data. 



6.3 Transients 

Transient AXPs are the most enigmatic objects. Most of them appear to be young sources whose 
field has not significantly decayed yet and their weak quiescent emission testifies of an extremely 
low efficiency in converting magnetic energy into X-rays. This could be related to magnetic dis- 
sipation (heat release) occurring only deep in the NS core, involving only a fraction of the whole 
volume and resulting in large neutrino energy losses and a low efficiency for the X-ray emission. 
As opposed to this, magnetic dissipation would be more distributed, or would occur closer to the 
NS surface, in persistent sources, thus reducing neutrino losses and leading to a higher efficiency 
of the X-ray emission. However, during outbursts transients become temporarily very similar to 
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the persistent sources and get much closer to them in the Lx - plane. The origin of this behaviour 
is still unclear and we will not discuss it further. 



7 DECAY OF THE INTERNAL MAGNETIC FIELD 



The presence of internal fields larger than the dipole components in magnetars has been consid- 
ered as a likely theoretical possibility since the suggestion made by TD96. In the previous sec- 
tions we have provided, for the first time, evidence based on the observed properties of persistent 
SGRs/AXPs that such a component must exist, if magnetic energy is indeed powering their X-ray 
emission. 

The decay of the internal field is even harder to constrain than the dipole, as the only obser- 
vational guidance is the evolution of Lx and the level of bursting activity (the latter being more 
qualitative in nature as it is harder to accurately quantify it). The X-ray luminosity depends on 
several physical details of the NS structure other than the properties of field decay. To keep the 
focus on the salient effect, we adopt a different approach here and, instead of dealing with general 
phenomenological decay models, we will calculate the expected evolution of Lx with Tc adopting 
two simple and general prescriptions from selected, physically-motivated models of field decay in 
NS interiors. 

There are two main possible locations for the decay of the internal fi eld. It could either take 



place in the liquid core of the NS, at 10^"^ < p < 10^^ g cm 



Thompson & Duncan 



2001 



Colpi et al. 



2000 



Arras et al 



rigid lattice of the inner crust, at p < 10^^ g crn ^ dVainshtein et al 



2001 



Arras et al. 



2004 



Pons & Geppert 



2007 : 



2004: 



(TD96, 



Hevl & Kulkami 



Dal 



Pons et al 



3). 



'Osso et al. 



2000 



1998 



2009h. or in the 



Konenkov & Geppert 



200% . In either case, heat released lo- 



cally by field dissipation is subsequently conducted to the surface, thereby powering the enhanced 
X-ray emisson. N ote that energy release may also take place at the NS surface, or just below it 



(p < IQio g cm 



Kaminker et al. 



2009L e.g. due to a gradual dissipation of electrical currents in 



a global/loca 



Beloborodov 



ized m agneto spheric twist (|Thompson et al 



2002 



Beloborodov & Thompson 



20071 : 



20091). This would likely allow a larger radiative efficiency. However, the twist and 



associated radiation would still draw their energy from that of an evolving, strongly twisted inter- 
nal field, either in the deep crust or core. 

As far as field decay in the liquid core is concerned, the state of matter there plays an important 
role. If it is normal npe matter, as opposed to superfluid, then ambipolar diffusion is expected to 
be the dominant channel through which the magnetic field decays. If, instead, either protons or 
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neutrons (or both species) were in a condensed state, par ticle interactions may be significantly af- 



fected, reducing or possi bly quenching the mode (GR92, 



Glampedakis et al 



Thompson & Dunca 



1996 



Jones 



20061 : 



201 ll) . We do not discuss the role of core condensation. Here we note that, even 



if ambipolar diffusion were completely quenched by core condensation, field decay in the crust, 
driven by the Hall effect, would still continue unaffected by the changing conditions in the core 



( |Arras et alj ([200^ give a quantitative account of this). 

Therefore, with the aim of illustrating just the salient effects on the X-ray luminosity of the 
decay of a strong internal field, we will consider here only the two limiting cases: ambipolar 
diffusion in a normal core of npe matter, or Hall-driven field decay in the inner crust, with a 
magnetically inactive core. The latter could be considered as a "minimal heating scenario " for 



magnetars. A study of re alistic mode 



dromagnetic instabilities (lArras et al. 



s, wh ich include the contribution to field dissipation of hy- 



2004 ) and magneto spheric currents ([Thompson et al. 



2002 



Beloborodov & Thompsonll2007h . along with the effects of strong crustal magnetic fields and light- 



element envelopes/atmospheres on radiative transfer, is clearly beyond the scope of this work. 



7.1 Field decay in the NS core 

We consider evolution of the internal fie 



to previous treatments (|Heyl & Kulkarni 



1998 



) through the solenoidal mode onl\ 



Colpi et al.l 



2001 



Arras et al. 



10 



. In addition 



20041), we allow ex- 



plicitly for a different decay law for the dipole field, according to our conclusions of sec.[5]and sec. 
[6l In our picture the decay of Sjnt heats the core and powers the surface X-ray emission, which 
will then decline following the decrease of the internal field. The decay of the dipole field will, on 
the other hand, determine the relation between real time, t, and spindown age, Tc (see Fig. [T]). We 
restrict attention to the two more realistic cases, a; = 1 or 1.5, as emerged in previous sections. 

For a core of normal npe matter the decay time for the solenoidal mode of ambipolar diffusion 
is Td^int ~ lO^Z-Sf/Ae y^^- The total available magnetic luminosity is then, according to Eq. |28l 
I^B,int = (i?,V3) Alt/rd,i„t ^ WB]^ll,Rl,eig s-\ 

As discussed in sec. |4l the equilibrium between heating and neutrino cooling determines the 
temperature as a function oi^ i?int, when ambipolar diffusion is active, as expressed by Eq. 
The evolution of the core temperature, Tc, will thus track directly that of fiint- 



The iiTotational mode evolves too slow to be of interest in the low-T regime and is very likely quenched by core condensation. 

We consider only modified URCA processes in the NS core. Note also that the equilibrium temperature in this case also depends on density, 

strictly speaking. However, we will carry out calculations at a fixed p = 7 X 10^* g cm~^, the average density of a 1.4 Mq NS with 10 Km 

radius, and focus only on the B-dependence. 
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Finally, an appropriate relation between the core and surface temperatures is needed in order 
to calculate the expected surface X-ray emission. This is the so-called Tb — relation, where Tb 
is the temperature at the base of the crust, the core being isothermal due to its large density and 
heat conduction. We adopt the minimal scaling for an unmagnetized, Fe envelope (Potekhin & 
Yakovlev 2001) 



% ~ 1.17 X 10*^ K 



5/4' 



1/4 



(29) 



where C = T^,9 - 0.001 9li\7 T^^gY^^ and ~ 1-87 is the surface gravity. 

Th e above relation is known to be sensitive to the strength and topolog y of crustal magnetic 



fields (iPage. Geppert & Ktikei 



2007 



Kaminker et al. 



20091 : 



Pons et al. 



20091) and also to the chem- 



ical composition of the outer layers of the crust. We comment later on these issues and their pos- 
sible relevance for our calculations. 

The surface luminosity, Lx = AuR^ashT^, is eventually obtained from Eq. |29l However, due 
to general relativistic corrections, an observer at infinity will measure the luminosity (Page et al. 
2007 and references therein) 
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which is the quantity we will compare with observations. The term in square parenthesis contains 
the dependence on i?int through the Tc-dependence of ^. 

Eq.|30] implies that fields in the 10^^ G range are strictly required to approach or even exceed 
~ 10^^ erg s^^, as observed in the youngest (Tc ^ 10^ yrs) sources. Note, however, that the 
total available power, L-q ^^^, would be 2-3 orders of magnitude larger during this early stage. 
Indeed, when magnetic energy is released in the core, most of it is carried away by neutrinos 
resulting in a very low efficiency of X-ray radiation. This can be estimated to be ex = Lx/LB,mt ~ 
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As the field decays, the equilibrium temperature in the core drops accordingly, the efficiency 
of neutrino emission decreases and ex grows. When it becomes close to unity, the NS thermal 
evolution becomes dominated by photon cooling and the equilibrium condition leading to Eq. [23] 
does not hold anymore. To treat this transition self-consistently we follow the evolution of 5int 
(and thus Tc and Lx) with time t from our chosen initial conditions. For a given choice of a and 
-Bdip.i, this is also known as a function of Tc. The temperature at which photon cooling becomes 
dominant, T,,, is defined by requiring that ex ^ 1/2, which is of course equivalent to finding 
where Lx first equals the neutrino luminosity. We denote by Lx,*, Td,*, Ant,* and t* the quantities 
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Figure 8. Maximum surface luminosity as a function of Tc of a NS with an internal magnetic field decaying through the solenoidal mode of 

o 6/5 

ambipolar diffusion, with ft: 9.6 X 10 Bj^^^ j^g yrs. The magenta arrows indicate the current upper limits on the spindown age and Lx of SGR 
0418+5729. Two initial values for the internal field are chosen, Bint.i = 4 X 10^^ G corresponds to red curves and 10^^ G corresponds to black 
curves. For each value of the internal field, three values of the initial dipole, B;jip,i = 10^^ , 3 X 10^* and lO'^* G were chosen, corresponding to the 
continuous, dashed and dot-dashed lines, respectively. Left Panel: The dipole field decays according to the o = 1 scaling, = 10^ yrs B~^.^ -^g; 
Right Panel ; The dipole field decays according to the a = 1.5 scaling, = 10^ yrs S^jp^j-. 

at this transition point. From here on, the heat released by the decay of the internal field will have 
to balance the energy lost to photon emission from the NS surface. That this equilibrium can be 
maintained is implied by the scaling Lx oc Tc while the total magnetic power is LB,int oc T'^- 
Hence, the latter decays (slightly) slower than the former as heat is lost and Tc decreases. 
The new equilibrium relation beyond thus reads 



with which we can eventually express the field decay time as a function of i?int 

rr«3.27xlO=yrsBrK?(^)"'" (32) 
showing that field decay becomes much faster now, with an effective a;int,ph = 2/19. 
Solutions for Bint, Tc and Lx as a function of real time are written straightforwardly 
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Lx,ooit-U) = . 20 (33) 
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and, for a given choice of a and i?dip,i, can be plotted versus the corresponding values of Tc. 

We stress that, in this regime, the former scaling LB,int oc 5;^/^ does not hold anymore. We 
obtained aint = 2/19 which implies LB,int oc Bf^/^ , matching the evolution of Lx as it must, given 
the equilibrium condition we imposed in the photon-cooling phase. 
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In Fig. [8] we show Lx.oo vs. curves for two different values of the initial internal field and 
three different values of the initial dipole. Two different choices for the decay index, a = 1 or 
1.5, are shown in the left and right panel of that figure, respectively. Note that curves for different 
values of -Bint,i and the same value of -Bdip,i become coincident, at late times. This happens because 
the surface luminosity tracks the instantaneous value of i?int and all models reach the same value 
of the internal field, at late times. On the other hand, curves with different i?dip,i and the same 
i?int,i maintain memory of the initial conditions although dipole fields also reach to the same value 
at late times. This happens because the relation between and t depends explictly on the initial 
dipole (cfr. Eq. [M]). At a given Tc, objects that had a different -Bdip,i are not co-eval. Those who 
had a weaker initial dipole are older, have a weaker Bi^^ and, thus, have a lower luminosity too. 

The relatively flat evolution of Lx for sources with < 10"^ yrs can be explained quite nat- 
urally, in this scenario. Larger initial fields produce larger luminosities (Eq. [30l) but, at the same 
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time, they have shorter decay times (rd,i oc -Bjnt ■ ). Hence, Lx is larger and bends earlier down- 
wards for increasingly strong initial fields, with the result that a flat, narrow strip of sources is 
produced. Its spread in Lx is smaller than the spread in -Bint,i and its maximum extension in 
roughly corresponds to the decay time, Td i of the minimum field. We estimate E^^^ff^ > lO^^ G. 

Two further properties of the above plots are worth noticing. First, the model with a faster 
decay of the dipole, a = 1, matches well the luminosities of persistent sources up to > 10^ 
yrs, but largely fails to account for the position of SGR 0418+5729. A direct link between the 
latter object and the persistent sources would then be impossible. Although the possibility that 
SGR 0418+5729 is linked to the transient sources is an interesting alternative, it has a serious 
drawback, if « = 1. A population of older (Tc ^10^ yrs), relatively bright magnetars would be 
expected, which is clearly not seen. All known sources at that (or beyond) are much dimmer 
and, accordingly, this possibility seems ruled out. 

The model with a; = 1 .5, on the other hand, provides a viable option to intepret SGR 041 8+5729 
as an older relative (tagc ~ 10® yrs) of the persistent sources whose dipole, as well as inter- 
nal, fields have strongly decayed. In particular, the two lower curves in the right panel of Fig. [8] 
give an internal field b[^^^^^ ^ 1.6(1.1), a dipole field -Bdip 1^2 ^ 6(4) and an X-ray luminosity 
-^x 3(f ^ ~ 30(0.4), respectively, from top to bottom. The predicted range of quiescent luminosities 
is quite wide, as opposed to the relatively narrow range for both Bi^t and -Bdip- This reflects the 
fact that SGR 0418+5729 is already in the photon-cooling dominated regime, where a steep drop 
in Lx occurs. A careful assessment of its actual quiescent emission would thus in principle put 
additonal constraints on the pysical parameters of this object. 
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7.2 Hall decay of the field in the inner crust 

The calculation of the previous section rules out small values of a and points tol.5<Q;<1.8as 
a viable option. Transition to core superfluidity (Tc < 5 x 10*^ K, Page et al. 201 1) is expected to 
occur at an age > a few x 10'' yrs, if i?int,i > 10^^ G, which for the viable values of a corresponds 
to Tf. ~ 10^ yrs. As stated before, we currently don't have a clear understanding of what will 
happen at this transition. The most conservative option is assuming that evolution of the core field 
will suddenly freez, its influence on the NS temperature soon becoming negligible. 

Even in this case, however, an internal field that threaded the NS crust would still be actively 
decaying, due to the Hall term in the indcution equation. We consider here this "minimal heating 
scenario" for magnetars, focusing on the effect of field decay in the NS crust. 

The timescale of Hall-driven decay of the magnetic field in NS crusts has a dependence on 
the actual field geometry, as was shown by Gumming et al. (2004). The field component that we 
are considering could either be a twisted (toroidal) field threading both core and crust (case /) or 
an azimuthal/multipolar field anchored only in the NS crust (case //). We consider the two cases 
separately. 



7.2.1 Case I Hall decay 

If the decaying crustal field were threading both the NS core and crust, its decay timescale would 
be sensitive to the global s tructure and is expecte d to be longer than Eq. |25]by a factor ~ R^/h, 
where h is crust t hickness iGumming et al.i (t2004,) . A more accurate expression was provided by 
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who also integrated the field induction equation through the crustal volume adopting this formula 
and a realistic density profile, to estimate the associated power output. The resulting expression 

LgL, ^ 2 X IQ^'Bl.R,,, erg s"^ , (35) 

exceeds the observed X-ray emission of younger sources if -Bint.i ^ 5 x 10^^ G, which accordingly 
represents a strict lower limit to the required crustal field. We evaluate the radiative efficiency of 
this model by considering, in a crude way, the impact of neutrino emitting processes within the NS 
crust. Neutrino bremsstrahlung and plasmon decay, in particular, become quickly very efficient in 
carrying away heat as the temperature rises, effectively limiting the maximum temperature that 
can be reached at the surface (cfr. TD96). An approximate, analytical expression for the implied 
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maximum surface te mperature, Ts™^'^^ , which also includes the effects of the magnetic field, was 



recently provided by (|Pons et al 



20091) 



^« - 1 + 0.002 LogBi2 y X iO K,^, ergs (35) 

where we neglect the very weak dependence on magnetic field in the last step. The effect on T^™'^^'^ 
would be much more pronounced if the field were predominantly tangential to the surface, thus 
strongly inhibiting heat conduction in the radial direction (Yakovlev et al. 2007; Page et al. 2007; 
Pons & Geppert 2009). However, the X-ray emission in the sample of persistent sources is very 
close to the limit of Eq.[36]and much higher than the limit derived for a strong tangential field. The 
effects of such a component do not appear to be relevant, at least for these sources, as was already 
pointed out by Kaminker et al. (2009). On the other hand, the effects of a strong tangential field 
would be qualitatively consistent with the weak quiescent emission of transient sources. A proper 
account of this issue is beyond our scope here and is postponed to future work. 

With the above formulae we can build approximate luminosity curves. As long as LB,int cal- 
culated through Eq. [35] is larger than Ls™^'''' (Eq. [361 ). our curves are limited by the latter value. 
Once -LB,int becomes comparable to L^^""^ the radiative efficiency is close to one and our curves 
track the evolution of LB,int from here on. The latter thus represents the bolometric luminosity of 
sources at later times. As in the previous section, we include the effect of gravitational redshift on 
the resulting Lx. Note that the Hall timescale is independent of temperature, so the evolution of 
-^B.int extrapolates from the previous, //-limited regime, into the photon-cooling regime. 

We calculated curves in the Lx,oo vs. plane in this way choosing two different values of the 
initial internal field, -Bint,i- For each value of Ant.i, two values of the initial dipole, -Bdip,i, were 
chosen, giving four curves in total for each value of the decay index a. Two different values of a 
were chosen and results for either choice are shown in the two panels of Fig. [9] 

Our main conclusions are summarized as follows: 

• the NS surface emission saturates at a level close to, but lower than, the observed luminosity 
of young AXPs/SGRs (Tc < 10^ yrs). Given the uncertainties and approximations implicit in the 
derivation of the limit in Eq.[36l the mismatch should not be regarded as a major issue. Further, as 
already noted, the gradual relaxation of a magneto spheric twist could provide an additional channel 
for release of the energy of the internal field directly at the NS surface. It is interesting to note that, 
despite not being limited by crustal neutrino processes in this case, a luminosity enhancement 
by just a factor of 2-3, at most, is all that is required to match the data, thus roughly confirming 
the overall energy budget estimated in our approximate calculations. A similar argument does not 
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Figure 9. Curves show the evolution of the NS surface luminosity, as a function of Tc , for Hall-driven decay of a field threading core and crust. 
The flat part of the curves con'esponds to the maximum surface luminosity of Eq. [36l and lasts as long as -LB.int is larger than this limit. Once it 
drops below that, curves switch to LB.int and, thus, represent the maximum available power from field decay {i.e. an upper limit to Lx.oo in the 
(2-10) keV range). Points mark the measured luminosities in the (2 — 10) keV range of magnetar candidates. The magenta arrows indicate the 
current upper limits on the spindown age and Lx of SGR 0418-H5729. Different symbols are explained in the figure. Left Panel: the dipole field 
decays according to the case a = 1 of sec.|5] with the scahng = W^B^^ yrs. Two values for the initial internal field are chosen, 3 X lO^*" G 
and 5 x 10^^ G. For each value of Bint.i, two curves are shown, conesponding to two different choices for the initial dipole field, 10^^ G (upper 
curves) and 10^** G (lower curves). Right Panel: same as left panel, but the dipole field is assumed to follow the o = 1.5 case of sec. [3] For each 
value of the internal field, three values of Sdip.i = 10^^ , 3 X 10^'* and 10^"* G are chosen. 



apply to the decaying parts of the curves, though, since those represent the evolution of L^^int, a 
limit that cannot be exceeded. 

• internal fields > 10^^ G are required also in this scenario. The minimal value i?int,i ^ 5 x 10^^ 
G provides a total magnetic power, L^-q\^^, marginally consistent with Lx,oo ~ 2 x lO^^reg of 
the youngest sources. Properly accounting for the powerful bursting activity of young sources and 
for realistic values of the radiative efficiency, would certainly imply a significantly larger minimal 
fiel(jl^. A neat example of this is provided by the 27 December 2004 Giant Flare fron SGR 1806- 
20, which released ~ 5 x 10^^ erg in high-energy photons. Even assuming that this was a unique 
event during the whole lifetime of the source, which we take to be Tc ^ 1.4 x 10^ yr, it would 
correspond to an average power output of ~ 1.5 x 10^^ erg s^^, an order of magnitude larger than 
the quiescent emission. 

• the internal field must also be able to provide its large power output for a sufficiently long 
time, as to match the duration of the apparent plateau in the -Lx,oo vs. up to > IC* yrs. This is 
comparable to the decay time of initial internal fields ~ 10^^ G. 

For a given value of i?int,i, a strong initial dipole and/or a small a push the end of the plateau to 
larger Tc, as discussed in sec. 17. II This is why the a = 1 model fails by overpredicting luminosities 

A self-consistent lower limit is obtained by setting the neutrino luminosity marginally equal to L^''"' . This implies LB,int twice as large and, 
flius, Bint,i > 7 X lO^'^G. 
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at Tc > 10^ yrs, unless dipole fields as low as 10^"^ G are assumed. Even for such a weak dipole, this 
model overpredicts the luminosity of SGR 0418+5729 and is thus ruled out. The a = 1.5 model 
seems, on the other hand, to well match the apparent bending at Tc ~ 10^ yrs, also remaining 
consistent with the position of SGR 0418+5729 (cfr. our estimate of its minimal power output in 
sec. [5]). 

Finally note that the value of Ts^^"^^ used in this section (Eq. [36l) would also apply in the case of 
core dissipation (sec. 17. II) . where we ignored this effect and let the surface luminosity free to grow. 
Despite this, the surface luminosities calculated in sec. 17. II are close to the maximum implied by 
Eq. |36]and very close, indeed, to the observed luminosities of magnetar candidates. In fact, also 
in that case the surface emission is limited by neutrino-emitting processes, which take place in the 
core rather than the crust. 

7.2.2 Case II Hall decay 

Finally, we consider the Hall-driven decay of a purely crustal field. In this case, the field would be 
completely insensitive to conditions in the core and its decay time is correctly given by Eq. [25l In 
a way completely analogous to Eq. [35l it is possible to define the total power output in the crust 
for this case as 

L^i;i ^ 2.7 X lO^'Bl.R,,, ergs-i . (37) 

Fig. [To] depicts the results the same calculation of the previous section, adopting the decay 
timescale Eq.|25]and the corresponding magnetic luminosity, Eq.[37l 

One conclusion can be drawn also in this case, which is in common to all other cases examined: 
an internal field of > 10^^ G is required to explain the observed luminosities of active magnetars, 
in particular middle-aged ones (tc > 10^ yrs). Indeed, although the luminosity of the youngest 
sources is orders of magnitude less than the available LB^mt at that age, luminosities in excess of 
10^^ erg s^^ are measured in middle-ages sources. This turns out to be a crucial property, since 
only fields larger than 10^^ G have a sufficiently large energy reservoir to account for it. 

For a = 1, it seems impossible to account for the luminosity of the AXPs 1E2259+586, even 
assuming the highest value of -Bdip,i that is consistent with the distribution of sources of Fig. [3] and 
m The AXP 4U 0142+60 is, on the other hand, just marginally consistent with this upper curve. 
Note, however, that both sources correspond to a lower initial dipole, for a = 1, as shown in Fig. 
[3l If they were both born with i?dip,i = 2 x 10^^, their spin periods would have to be < lis, as 
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Figure 10. Curves show the evolution of the NS luminosity as a function of Tc, for Hall-driven decay of a purely crustal field. The flat part of 
the curves corresponds to the maximum surface luminosity of Eq. [36] and lasts as long as Lb int is larger than this Hmit. Once it drops below 
that, curves switch to Lb, int and, thus, represent the maximum available power from field decay {i.e. an upper limit to Lx.oo in the (2-10) keV 
range). Points mark the measured luminosities in the (2 — 10) keV range of magnetar candidates. Different symbols are explained in the figure. The 
magenta arrows indicate the current upper limits to the spindown age and Lx of SGR 0418-1-5729. Left Panel: the dipole field decays according to 
the case a = 1 of sec. [5] with the scaling = 10*^ yrs. Two values for the initial internal field are chosen, 3 X lO"' G and 10^'' G. For each 
value of the internal field, three curves are shown, coiTesponding to three different values of the initial dipole field, 2 X 10^^ G, 5 X 10^^ G and 
2 X 10" G. Right Panel: same as left panel, but the dipole field is assumed to follow the o = 1.5 case of sec.[5]and different representative values 
of i3dip,i chosen. 



clearly seen in Fig.lH For -Bdip,i ^ 10^^ G the mismatch of model curves with observations would 
be far too large. This discrepancy is thus very significant. 

A wide range of quiescent luminosities is predicted for SGR 0418+5729, essentially depending 
on its initial dipole field. The latter is likely close to -Bdip,i ~ 10^^ G for a = 1, (see Fig. [3] and ID), 
for which a quiescent luminosity ~ 2 x 10'^^ erg s"^ would be expected. A direct measurement of 
its quiescent emission may therefore provide a conclusive answer about this scenario. 

The case a = 1 would have to be considered at best marginal, given the above considerations, 
essentially for the same reasons as the case « = 1 in sec.[6l The additional general arguments we 
presented against this small a value make this scenario, in conclusion, very unlikely. 

The case a = 1.5 falls short of the X-ray emission of several sources beyond Tc > 10^ yr and 
is, therefore, completely ruled out. 

Overall, unless a different scaling for rnaii can be provided, or unless an additional mechanism 
for decay of purely crustal fields is proposed, this hypothesis appears inconsistent with the obser- 
vational properties of the source sample considered. Accordingly, we are led to rule it out in its 
present form. 
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8 SUMMARY AND CONCLUSIONS 

In this work we have studied the implications of dipole field decay on the dynamical evolution of 
a spinning, magnetized NS, assuming a power-law scaling for the field decay time, oc B'^^ and 
i?dip oc -B^J". If the field decays sufficiently fast (a < 2) the spindown time grows faster than 
the decay time, r^. Once > r^, the dipole continues to decay while the NS spin period hardly 
changes and an asymptotic value of the spin period. Poo, is reached. For a slow (a > 2) decay, the 
above condition is never met and the NS continues spinning down forever. 

The distribution of magnetar candidates (SGRs, AXPs, transients) and XDINs in the B^ip vs. 
Tc ( or Bciip vs. P) plane provides a strong evidence for a fast {a < 2) dipole field decay. This 
conclusion is based on two striking properties. First, the absence of objects with large periods, 
well above 10 s, and in particular with large dipole fields B^ip and large spindown ages Tc. We 
know that objects with large B^ip exist since we find them at small Tc. Similarly we know that 
objects with large Tc also exist, but they are all found only with relatively small -Bdip- If -Bdip did 
not decay then = t and there would be no reason why there are no old (large r^) objects with 
large -Bdip- Note that such objects should be more numerous that younger ones and they should be 
at least as detectable as objects with a similar and smaller -Bdip- 

Decay of the dipole field offers the most natural interpretation of this fact, as it naturally pro- 
duces bending evolutionary tracks in the i?dip vs. Tc plane or, equivalently, the P vs. P and Pdip vs. 
P planes (see Figs.[3l|4]and[5]). In particular, if the decay law has a < 2, an asymptotic spin period 
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is approached which corresponds to the scaling Pdip,i5 ~ Poo.io kyr • Sources that reach their 
asymptotic period move along a trajectory with a constant period. Poo, which, in turn, is related to 
the value of a and of the initial magnetic dipole, Pdip,i (see Eq. fTT]) . 

The second key feature is the fact that spin period periods of all (but one) sources with ^ 
10"^ yr are within a very narrow strip corresponding to 7 < Poo < 12 s. This supports the existence 
of an asymptotic period and hence of a single decay mechanism operating in all sources, with 
a < 2 and with a decay time ~ lO'^ Pfg" yr (sec. [5]). The narrowness of the distribution in Poo 
implies that the decay index must be in the range 1 < a < 2. 

Hall decay in the crust, corresponding to a = 1, provides a reasonable explanation for the 
observed Pdip-Tc distribution. It implies that all SGRs, AXPs, transients and XDINs were drawn 
from a distribution of Pdip,i G [0.2 — 2] x 10^^ G. A similarly good description of the data is 
provided by larger values of a, for which we consider a = 1.5 as representative. These allow for 
a wider distributions of initial dipoles, making them favored in this respect. As a approaches 2, 
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asymptotic periods > 12 s would necessarily be expected and the current data disfavor «-values 
too close to 2. Note that all viable models imply that there cannot be a population of sources whose 
dipole field at birth is significantly larger than 2 x 10^^ G. These would indeed reach asymptotic 
spin periods significantly longer than observed. 

An important implication of the dipole field decay with a < 2 is that the sources are younger 
than what their spindown age implies {t < Tc). The current value of t/r^ also depends on i?dip,i- 
For magnetar candidates this ratio is significantly smaller than unity and their real ages can be a 
factor of ~5-50 smaller than t^, depending on their initial field, age, and exact value of a. 

The ratio I/tc decreases for decreasing values of a, at a given Tc (see Fig.[T]). For a fixed value 
of a,t/Tc decreases as the object ages. 

X-ray luminosities of XDINs, their clustering at Tc ~ a few x 10^ yr and the distribution of their 
effective temperatures point to their ages being t > 10^ yr. This further constrains possible values 
of a. Specifically, a = I decay implies that these sources are younger, or even much younger, than 
this value. On the other hand, values of a ~ 1.5 well match this age requirement for XDINs. Thus, 
the case a = 1 is basically ruled out and a more restrictive condition, 1.5 < a < 1.8, appears to 
be the most consistent with the whole available data. 

Although the decaying dipole field dissipates energy, it turns out that unlike the simplest ex- 
pectations, this energy is insufficient to power the strong X-ray emission of SGRs/AXPs (and, 
possibly, other related classes). While exponential decay is rapid enough that L^^dip would suf- 
fice for all sources, the current upper limit on the weak quiescent emission of SGR 0418-1-5729 
(< 6 X lO'^^ erg s~^) already implies « > 1, as otherwise this source would be too young to be so 
dim. 

For a = 1, the maximal available power, Le.dip, falls short of the observed X-ray emission 
of persistent AXPs at Tc ~ 10^ yr and the current observed emission of SGR 0418-1-5729 is far 
largeiLjthan its i^B.dip for all values of a ^ 1. Hence, a = 1 and energy provided just by magnetic 
dipole field decay should be considered at best marginally consistent with the observed distribution 
of sources in the Lx vs. Tc plane. However, as we have discussed above, a = 1 is ruled out by the 
requirements on the ages of XDINs, which altogether imply 1-5 < a < 1.8. For such values of a 
(or even for a > 1) it is impossible to power the X-ray emission of SGRs/AXPs by the decay of the 
dipole field. If magnetic energy powers this emission, the available reservoir must be larger than 



This post burst luminosity might over estimate the real quiescent emission. However even the emission from the outbursts of SGR 0418+5729 
suggests a minimal energy output larger than Lb, dip. for this source. 
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that of the dipole and must decay on a somewhat longer timescale. The most natural conclusion is 
that this is provided by the magnetic field in the NS interior (sec. (6]). 

In order to identify at least the basic properties of the internal field, we considered its decay 
either in the liquid core, through the solenoidal mode of ambipolar diffusion, or the inner crust, 
through the Hall term. In the latter case, we treat separately fields threading the whole NS volume 
or fields confined to the crust alone. A general conclusion in all cases is that the observed large Lx 
values at up to Tc ~ 10^ yrs require a very large initial internal field, -Bint,i ^ 10^^ G . Combined 
with the maximum value of ~ 2 x 10^^ G for the initial dipole field this implies that, at formation, 
the internal magnetic field in magnetars contains at least 20-30 times more energy than the dipole, 
likely one order of magnitude more than that. This raises a question about the stability of the 
implied configurations, which is a currently open field. Although it is difficult, at this stage, to 
have a precise answer on global field stability, recent investigations (cfr. ?) have highlighted the 
key role of stable stratification of stellar interiors in providing stability to predominantly toroidal 
fields. As opposed to this, it was shown that a predominantly poloidal field would remain unstable, 
independent of the interior stratification. The maximum ratio of toroidal-to-poloidal field strength 
that can be stablized is estimated, conservatively, to be < 10, for a ~ 10^^ G toroidal field, roughly 
consistent with our results. Clearly, more work is needed to settle this important topic. 

Fields threading the NS core have relatively long decay timescales (> 10"^ yr for B-^^t = 10^^ 
G), whether their dissipation occurs primarily in the crust or the liquid core. This implies that 
it takes > 10"' yr for Lx to decrease appreciably. On the other hand, observations show that this 
decrease occurs at Tc < 10^ yrs, approximately 10 times larger than the real age. This favors values 
of a > 1 and, in particular, our choice a = 1.5 meets easily this requirement (cfr. Fig.[T]). If a ~ 1, 
on the other hand, Tc grows too fast with real time and when Lx eventually starts decreasing, Tc is 
already too large, in sharp contrast with observations. 

If the internal field is confined to the NS crust, on the other hand, it decays on a ~ 10^ yr 
timescale and the situation is reversed. The case a = 1.5 largely fails because it implies that a 
given Tc is reached when the source is too old and its X-ray emission too weak. The case a = 1 is 
marginally consistent with Tc ~ 10^ yr-old sources, and only for an extreme choice of parameter 
values. Considering the strong independent arguments against such a low a value, we conclude 
that purely crustal fields cannot account for the distribution of SGRs/AXPs in the Lx vs Tc plane. 

SGR 0418-1-5729 has the largest Tc, weakest B^ip and Lx of our sample and thus provides 
crucial insight into the long-term evolution of magnetar candidates. It is thus of particular interest 
to summarize our main conclusions on its properties. For the reasons just mentioned, we restrict 
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attention to the two scenarios where the internal field threads the whole NS volume and, thus, to 
dipole field decay with a ^ 1.5. We find very similar conclusions in the two cases, whether the 
internal field decays in the crust or liquid core. This is not surprising, given the very similar values 
of Oint and ra^int in both models. From the decay of the dipole field with a = 1.5, this source was 
likely born with -Bdip,i ~ (3 — 5) x lO^'' G, which implies its real age is tage ~ (1 ~ 2) x 10^ yrs, the 
youngest age corresponding to the larger initial field. With these assumption on the initial dipole, 
its current dipole field would be ~ (4 - 7) x 10^^ G and the internal field Bi^t ~(l.l-2)xl0^'^ 
G, independent of its initial value (as long as it was > lO^^G). The average luminosity due to the 
decay of the internal field is LB,int ~ Lx ~ (4 — 10) x 10^° erg s"^. 

Note that the internal field is expected to be close to the lower limit for being able to produce 
crust breaking (and, thus, bursts, according to the current understanding of these events, cfr. TD96; 
Pema & Pons 2011). Also note that the expected luminosity can be checked by continued mon- 
itoring of this source and is, incidentally, similar to the estimated average output from outbursts 
(cfr. sec. [5]). 



9 DISCUSSION & EVOLUTIONARY SEQUENCE 

In the light of our results, we can sketch a tentative evolutionary picture for the different classes of 
o bjects considered here. The most obvious link is between SGRs and AXPs, as already suggested 



by lKouveliotou et al. 



(119981) . 

These sources appear to form a continuous sequence, with persistent SGRs seen as the youngest, 
mostly magnetized and brightest sources. Persistent AXPs are found at comparable field strength 
and similar, although slightly older, ages. A couple of evolved AXPs are known (1E2259+586 
and 4U 0142+61), whose dipole fields have decayed substantially and whose X-ray luminosity is 
somewhat weaker than in younger objects. The X-ray luminosity of SGRs/ AXPs is inconsistent 
with the decay of the dipole field alone and needs to be powered by an additional, stronger field 
component. The latter must thread the whole NS volume and its strength likely exceeds 10^^ G at 
birth. At young ages the internal field likely decays through a combination of ambipolar diffusion 
in the core and Hall decay in the crust. The released heat powers the strong thermal emission typ- 
ical of these sources. Neutrino-emitting processes, whether in the liquid core or in the crust, limit 
the thermal luminosity of young magnetars to be < 10^^ erg s^^ and provide a likely explanation 
for the observed plateau in the distribution of Lx up to Tc > 10"^ yrs. 

Transient SGRs/AXPs are not related in an obvious way to persistent sources. Their dipole 
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fields are systematically weakeo and their spin periods also appear to be shorter, on average, 
suggesting that most of them have not yet reached the asymptotic spin period. The two oldest tran- 
sients have significantly weaker dipole fields and somewhat slower spins (CXO J164710-455216 
and SGR 0418+5729), thus supporting this basic picture. 

The decay of Sdip appears consistent with having the same properties as in persistent sources. 
However, the X-ray luminosities of transients are much weaker and it is not possible, at present, 
to identify a decline at ~ 100 kyrs, as seen in persistent sources. 

The weak quiescent emission of transients is very striking, being much lower even than the 
power provided by their decaying dipole field. The relatively bright, short-time bursts they emit in 
outburst and the ensuing large increase in X-ray luminosities, that make them temporarily similar 
to the persistent sources, supports, on the other hand, their connection to the latter groujolf this is 
indeed the case, then, it is possible that transients also host a decaying, strong internal field, despite 
the very low radiative efficiency of their quiescent state. We can only speculate, at this stage, that 
a strong azimuthal component of the crustal field provides a possible explanation for their X-ray 
deficiency. If the energy released by magnetic dissipation is confined to a sufficient depth in the 
NS interior, thermal insulation of the hot regions by such a field component could suppress the 
surface X-ray emission by a large factor (cfr. Potekhin & Yakovlev 2001, Kaminker et al. 2009). 
The difference with persistent sources would then be attributed to the existence of a much more 
ordered field component in transients, which globally suppresses the radiative efficiency ex. Only 
at outbursts does the internal field gain access to the more efficient channel of X-ray radiation, 
probably due to flaring of smaller-scale field structures temporarily creating a preferential path 
for heat conduction, or a localized dissipation region close to the NS surface. This is qualitatively 
consistent with the small (and shrinking) emitting areas and high (and decreasing) temperatures of 
the extra BB spectral components found in outburst emission from transients. 

In this context, the dichotomy in behaviour between transient and persistent SGRs/AXPs which 
likely represent two distinctevolutionary sequences (see Fig. fTT)) . could point to a different mecha- 
nism for generation of their strong fields at formation. One might speculate, based on the apparent 
properties of their magnetic fields, that SGR/AXPs appear to be bom with somewhat higher ini- 
tial dipole fields (-Bdip.i ~ 10^^ G) compared to transients (which are narrowly clustered around 
5dip,i ~ 2 X 10^^ G). SGR/AXPs also have a strong initial internal chaotic field (5int,i > 10^^ G) 

^'^ It is remai'kable that most of them have -Bdip ~ 2 X lO^** G. 

During outbursts, the X-ray luminosity of transients largely exceed their rotational energy losses, like in persistent sources, which is not the 

case in quiescence (cfr. Tab.[T). 
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that varies significantly on length scales much smaller than the NS radius, which naturally ac- 
counts for their good heat conduction from the core to the surface and for the more frequent and 
more powerful bursting activity. Transients, on the other hand, likely have a strong and globally 
ordered tangential field, which efficiently suppresses heat conduction from the core to the surface 
and also results in much fewer and smaller bursts (that are likely powered by small scale field 
structures). Thus, one might speculate that SGR/AXPs might attain their fields at birth through the 
alpha-omega dynamo, if they were initially sufficiently rapidly rotating, while NSs in the transient 
branch might obtain their fields through a different (unclear as of yet) mechanism. 

The -Bdip vs. P distribution of XDINs is also consistent with the same dipole field decay of 
other classes. We conclude that these sources represent ~ (1 — 6) x 10^ year-old NSs bom with 
5dip,i ~ (0.03-2) X 10^^ G (for a = 1.5; see bottomleft panel of Fig.©. Theirinitial dipole fields 
have significantly decayed, making them apparently look older than they really are. Hence, their 
weak X-ray emission is likely dominated by their remnant heat, although a minor contribution from 
the decay of their dipole field cannot be ruled out, in general. As opposed to other classes, they 
show no evidence for an additional, stronger field component in their interior. As such, their most 
likely connection is with other high-field NSs not showing any peculiar magnetar-like activity. The 
brighter XDINs had initial dipole fields comparable to those of SGR/AXPs or transients, but they 
likely have a weaker internal field. The dimmest and fastest spinning of all XDINs, RX J0420, 
was bom with the weakest dipole field, -Bdip,i ~ 4 x 10^^ G (for a = 1.5). This is well below that 
of SGR/AXPs, or even of transients, and actually suggests a possible evolutionary link with either 
the high B-field radio pulsars or the Rotation-Powered X-ray Pulsars (RPPs), which would depend 
on the strength (and, possibly, the geometry) of both the initial dipole field and the initial internal 
field. This tentative optiortjis also sketched in Fig. [TH alongside the clearer evolutionary tracks 
of the SGR/AXP branch and the transient branch, that have been discussed above. 

How does the "weak field magnetar" SGR 0418-1-5729 fit in the general picture? This object 
was likely bom with a dipole field ~ (3 — 5) x 10^^ G that has now decayed to (4 — 7) x 10^^ G, 
and an internal field > 10^^ G, which has now decayed to ~ (1 — 2) x 10^^ G. Despite its current 
appearance as a transient source, SGR 0418-1-5729 has most likely followed the track of persistent 
SGRs/AXPs in the Lx vs. Tc plane. The absence of more luminous sources at similar Tc implies that 
the emission of persistent sources must drop significantly in the range Tc ~ (10^ — 10'') yrs. It is not 
clear whether the decrease in luminosity of the persistent sources passes through the upper tail of 



Data for the RPPs are taken from Tab. 3 in Mereghetti (2011). 
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Figure 11. A schematic diagram of the evolutionary sequences we infer for different high magnetic field NS source classes in the -Bjip - f"c plane 

\f'2 ' 1/2 

(which can be thought of as a rotated version of the P-P plane since P oc -Bdip^c and P Sdip''c )■ Members of different classes are 
denoted by different symbols, while the classes as a whole are indicated by colored ellipses. We clearly identify two distinct evolutionary tracks - 
the SGR/AXP branch (solid arrow) and the transient branch (short-dashed arrow), which start with different initial dipole fields (and likely also 
different internal fields) but converge to similar values (at a given Tc) after their fields decay significantly. The dimmest and fastest spinning XDIN 
(RX J0420) must have had a significantly smaller initial dipole field than either of these two branches, and thus must have followed a different 
evolutionary track, which might be related to rotation-powered X-ray pulsars (RPP) or high-field radio pulsars (this tentative track is indicated by 
the dotted arrow). Such a significantly lower initial dipole field results in a somewhat faster asymptotic spin period. 

XDINs (like the lower curves in the right panel of Fig. |9]may suggest) or completely avoids them 
(like most curves in the right panels of Figs. [8] and |9l). In the former case XDINs would be a less 
homogeneous class of sources, which included both passively cooling, high-field NSs and some 
old magnetars. This would also imply that SGR 0418+5729 follows the brightest XDINs along the 
same cooling sequence. Its quiescent luminosity should thus be orders of magnitude lower than 
that of the brightest XDINs. In the latter case, the XDINs are just passively cooling NSs, with no 
particular relation to persistent SGRs/AXPs, accidentally occupying a region in parameter space 
that is between persistent and transient magnetars. 

Overall we have demonstrated that the dipole magnetic field of magnetars must decay on a 
time scale of (10^ /B^-^^ years. The most likely decay law has 1.5 < a < 1.8. At the same time 
we have shown that the power supplied by the decaying dipole field is not sufficient to power the 
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X-ra y luminosity of these objects. Un 



(cfr. 



Triimper et al. 



2010 



Alpar et al. 



ess th ere is an external energy reservoir, like a fallback disk 



20111), this energy can arise only from a stronger internal 
field whose decay is decoupled from the decay of the dipole^. Detailed analyses of the X-ray 
emission depends on the structure of this internal field, its decay mode as well as cooling and heat 
transfer within the NS. Further observations of this exciting group of object could reveal some of 
the hidden secrets of the NS interiors as well as provide clues on magnetic field generation and 
amplification at their birth. 
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